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Suinnary 

This report describes NBS work for NASA In support of NASA's Materials 
Processing In Space Program under NASA Government Order H-27934B (Proper".les 
of Electronic Materials) covering che period April 1, 1980 to March 31, 1981. 

The work emphasizes measurement of materials properties and thermophysical 
properties Important In materials processing applications. It has two main 
thrusts: 

1) Carrying out precision measurements In space and Investigating the 
feasibility of Improved measurements when the space environment offers a 
unique opportunity for performing such measurements. These measurements 
would be useful for either space processing or processes on the ground. 

2) Obtaining precision measurements on materials properties when these 
properties are important to the design and Interpretation of space processing 
experiments. These measurements would be carried out either in space or on 
the ground. 

This work has been carried out in five tasks. These tasks have, as two 
of their focal points, the role of convection effects and the role of con- 
tainer effects, both of which would differ in space-based experiments from 
those found in ground-based experiments. The results obtained for each task 
are given in detailed summaries in the body of the report. 

Briefly, in Task 1 - Surface Tensions and Their Variations with Temperature 
and Impurities - special equipment has oeen put into operation which, for the 
first time, allows simultaneous liquid-vapor surface tension and Auger surface 
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concentration measurements to be made on metal liquids* This apparatus has 
been used to analyse surface tensions and reactions on pure gallium and a 
gallium-tin alloy » providing information which is important for prediction 
and analysis of surface processes on liquids in space flight experiments* 
Applications include, for example, Marangonl convection and surface segregation 
effects occurring during solidification* In this work, ionic bombardment 
sputtering with argon was found to be a very promising method for producing 
clean liquid surfaces, since unexpected convection currents appeared which 
allow rapid and efficient cleaning of carbon particles from the wuole surface, 
or indeed the whole volume, of the liquid droplet* For gallium specimens 
cleaned by this method, a linear dependence of surface tension on temperature 
was found, in contrast to the quadratic dependence found for less clean 
samples* For Ga-Sn (1*9 at*%) alloys the Auger measurements showed strong 
temperature-dependent surface segregation of Sn to the surface* This result 
was consistent with the simultaneous surface tension measurements made on the 
alloys* 

In Task 2 - Convection During Unidirectional Solidification - the 
succinonitrile-ethanol system is being investigated as a transparent material 
in which convective flows produced by simultaneous temperature and composition 
gradients during unidirectional solidification can be quantitatively 
measured and compared with theoretical predictions* Knowledge of convective 
conditions is Important in comparing solidification in microgravity with that 
on earth and in evaluating the advantages of space processing* In order to 
determine the distribution coefficient k, the phase diagram in the low 



ethanol range was determined, and k was measured as 0.044. Calculations to 
predict the resume ot convective Instabilities then were made using this 
distribution ci ef f icient value. It was predicted that ethanol concentrations 
in the range 10 to 10 weight percent can provide either convective, 
inter facial or oscillatory instability depending on the values of temperature 
gradients and solidification velocities. An apparatus has been constructed 
to measure convective flows in these ranges in this material, and initial 
baseline thermal convection measurements have been made. Velocities are 
determined by time lapse photography of polyvinyltoluene spheres having 
nearly the same density as the solidifying liquid with the spheres acting as 
markers in the liquid to map the convective flows. 

In Task 3 - Measurement of High Temperature Thermophysical Properties of 
Tungsten Liquid and Solid - this task contributes to a joint effort of the 
General Electric Advanced Applications Laboratory (GE), Rice University 
(RICE) and the National Bureau of Standards (NBS) to develop levitation/drop 
calorimetry techniques. The immediate objective is to measure the specific 
heat of liquid tungsten, in order to provide key benchmark measurements on 
thermophysical properties. The current phase of thi.s effort has involved 
completing the physical connection of the RICE calorimetric system to the GE 
levxtation facility. In this phase, improvements have been made particularly 
in the temperature measurement area* A data logging system has been added to 
eliminate the need for a human obaerver of the calorimeter heat trend. This 
automation, along with proposed automation of the calorimeter data system, 
which is made necessary by the erratic levitation holding times of the molten 
tungsten drops, can also be expected to aid in automatic operation of similar 
space flight experiments. 
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In Task 4 - High Temperature Thermophysical Properties of Refractory 
Materials from Free Cooling Experiments - a free-cooling method for measuring 
specific heats of reactive and refractory materials is being developed. This 
method allows use of the pcsitionlng-levitation capabilities available in 
space to provide thermophysical measurements without contact with a container. 
Tests with a millisecond pulse heating system have been performed using pyro- 
electric radiometer de^’ices to measure the energy radiated from the specimen. 

In addition, calculations have been performed to relate the surface temperature 
and radiated energy of spherical samples (which can be measured directly with 
a pyrometer and the newly developed pyroelectric radiometer devices) to the 
specific heat, thermal conductivity, and emissivity of the sample. Under 
certain conditions, it should be possible to determine the temperature depen- 
dence of these parameters from simultaneous measurements of energy radiated 
and temperature change during free-cooling of a levitated sample. 

In Task 5 - Experimental and Theoretical Studies in Wetting and Multilayer 
Adsorption - elllpsometric parameters of the liquid-vapor interface between a 
binary mixture and air have been measured near the predicted 

trsi^sition from high adsorption to low adsorption. Also a method of measuring 
and recording the index of refraction of this transparent binary mixture 
simultaneously with measurement of the elllpsometric parameters has been 
tested. Separate laser light sources are used for the two types of measure- 
ments. Since the index of refraction measurement is e>:fremely sensitive to 
the appearance of macroscopic amounts of a second liquid phase, it can serve 
to precisely locate the miscibility gap. Thus the simultaneous measurement of 
both the index of refraction and the elllpsometric parameters will locate the 
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adsorption transition with respect to the miscibility gap in the same sample* 
Measurements of ellipsometric and other optical parameters will allow evaluation 
of wetting behavior involving immiscible materials, a topic of major significance 
in space processing applications. 
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Task 1 


Surface Tensions and Their Variations 
with Temperature and Impurities 

S. C. Hardy 
Metallurgy Division 
Center for Materials Science 

and 

J. Fine 

Surface Science Division 
Center for Chemical Physics 

Summary 

Sessile drops of pure gallium and a gallium- tin (1-9 at.%) alloy have 
been studied in an Auger spectrometer. The surfaces of the drops were cleaned 
by heating to 500 ®C and by ionic bombardment with argon- The surface tension 
of pure galliuHi cleaned in this way was found to vary linearly with tempera- 
ture and could be represented by 

Y * 724.1 - .0721 T (®C) , 

This linear relationship is significantly different from the quadratic 
variation found previously for samples cleaned only by heating in vacuum and 
hydrogen. The difference in the two results is tentatively attributed to the 

presence of small carbon platelets on the surfaces of the unsputtered drops. 

« 

The segregation of tin at the surface of gallium-1. 9% tin drop was 
measured by comparing the tin Auger line amplitudes with that of a pure tin 
sample with the same spectrometer settings. Although there is some uncertainty 
in the attenuation length and backscattering factors used in the analy:*’is. 


Preceding page blank 
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surface concentrations in excess of two monolayers were found at 30 ®C, w'tb a 
sharp decline with increasing temperature. Surface concentrations calculated 
from surface tension measurements were in general agreement with this conclusion. 
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Introduction 


It has been widely recognized that surface tension measurements of 
liquid metals and semiconductors are flawed by the absence of any independent 
characterization of the chemical species present on the surfaces. The devel- 
opment of medem surface spectroscopies makes available a host of techniques 
which can address this gap in our knowledge of liquid surfaces, at least for 
some low vapor pressure materials. In previously reported work. Auger electron 
spectroscopy (A.E.S.) was applied to liquid gallium and provided valuable 
information about t le state of the surface. The configuration of the spec- 
trometer, however, did not allow the surface tension to be measured with the 
desired accuracy. This past year has been largely spent in the assembly and 
operation of an Auger instrument which does permit accurate surface tension 
measurements of the liquid sample. 

Our initial studies of liquid gallium last year showed that the primary 
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surface Impurities were oxygen and carbor^ ’ * . The oxygen was undoubtedly 

in the form of a gallium oxide and could be effectively redUwCd and ranoved by 
heating; the carbon, however, remained in the form of small platelets floating 
on the surface. In an attempt to clean a spot on the surface to obtain the 
spectrum of clean gallium, we sputtered a small area with argon ions using a 
conventional low power gun. The ion bombardment was found to generate large 
scale fluid flow that eventually carried all of the carbon into the ion beam 
where it was sputtered away. The result was a gallium surface which appeared 
to be perfectly clean within the measurement limits of the Auger spectrometer • 
Our first surface tension measurements on gallium sessile drops v7ith surfaces 
cleaned by sputtering will be seen to have a different temperature dependence 


frcm our previous results with thena^lly cleaned surfaces* Tills observation 
reinforces our feeling that the sp i er cleaning process will be an important 
tool in future studies of liquid surfaces. We will present here our current 
ideas on the meclianism responaiblo for the fluid flow uuring the sputtering 
process • 

The study of surface segregation in liquid alloys is an area which seems 
partlcul .jly suited to A.E.S. The comparison of directly measured surface 
concentrations with those deduced from surface tension values should eventually 
result in improved models for the interface. We have begun measurements with 
dilute alloys of Sn in Ga and will describe some Initial results. 

Experimental Apparatus and Procedures 

The Auger spectrometer which we have assembled is a conventiona?. single- 
pass cylindrical mirror analyzer system with a coaxial electron gun. The 
CMA is mounted vertically so that it can be focussed on the apex of a 
sessile drop located beneath it. This arrangement was chosen because it 
maximizes the fraction of the drop surface which can be analyzed by the 
spectrometer and provides excellent visibility for photography by locating 
the drop near the center oi two opposing 8" view ports. The shallow cylin- 
drical cup containing the sessile drop is clamped to a simple hot stage 
mounted on a horizontal manipulator which translates in three orthogonal 
directions. An ion gun mounted in a side port is also focussed at the drop 
apex. Fig. 1 shows schematically the experimental arrangement. The stain- 
less steel bell jar in \Aiich these components are mounted is equipped w'th 
three pneumatic vibration isolation supports which are adjusted to level 
the cup for sessile drop surface tension measurements. The bell jar is 
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evacuated with sorption and ion pumps to minimize hydrocarbon concentrations* 
Pressures in the mid 10 torr range are achieved without baking as measured 
by a nude ionization gauge. High purity gases can be admitted to the bell 
jar from a connecting system for sputtering or oxidation studies* 

The hot stage is also shown schematically in Fig* 1* The seven heaters 
are 0*6 mil tungsten wire spirals enclosed in alumina tubes and are pressed 
to the bottom of the plate which supports the cup by an opposing plate* The 
heater ends are shielded to prevent evaporative contamination of the sessile 
drop and the CMA. The cups are generally made of 0*6” O.D* silica cylinders 
fused to 1” diameter flat silica plates to provide a wide shoulder at the 
base. The body of the cup passes through a hole in a thin steel sheet which 
presses the cup base against the heated plate* A chromel alumel thermocouple 
is mounted on this pressure plate at the edge of the hole through which the 
cup passes. The temperatures read by this thermocouple were in agreement 
with those measured by a W-Rh thermocouple inserted into a sessile drop of 
gallium contained in the cup. 

The hot stage, which is constructed of stainless steel, is mounted by 
three thin screws to a silica plate* This silica plate is itself attached to 
the platform of the manipulator and provides effective thermal isolation for 
the hot stage by its length and its low conductivity. Consequently, the 
equilibrium temperature of the stage is extremely stable* This thermal 
isolation, however, imposes the penalty of a slow cooling rate at low 
temperatures where radiation losses are small. Fig. 2 shows the equilibrium 
temperature of the hot stage as a function of heater current. The power 
source was a small variable transformer. 
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The use of insulating cups to contain the sessile drop requires that an 


electrical contact be made with the sample for Auger analysis. If this is 
not done, charge builds up on the sample and disrupts the operation of the 
CMA. It is also Important to be able to measure the electron beam current in 
order to compare different Auger spectra. Wc make electrical contact to the 
drop with a fine W wire spot welded to a nickel rod that is rigidly held in 
an electrical feedthrough on the manipulator flange. Ft-r surface tension 
measurements, the manipulator is used to translate the drop and break contact. 
This is necessary to obtain an unperturbed sessile drop. This technique was 
also used to Insert a fine thermocouple in the drop as mentioned previously. 

Prior to forming a sessile drop, the bell jar was evacuated and the hot 
stage temperature was raised to 500 "C for about an hour. After cooling, the 
system was back filled with dry nitrogen gas. A small port was removed and 
gas allowed to flow through the system while the cup was partially filled 
with liquid gallium using a syringe. The port was then replaced and the 
system very slowly evacuated with the sorption pumps. If evacuation is too 
rapid the expansion of gases trapped under the gallium in the cup will expel 
the liquid charge. After evacuation the process is repeated to form a 
sessile drop large enough for surface tension measurements. Fig. 3 is a 
photograph of a sessile drop formed in this way. 

The instrumentation for our Auger spectrometer is not yet complete. The 
electronics used to form a video image of the sample by rasterlng the electron 
beam over the surface were not received and have had to be reordered. The 
imaging of the sessile drop is important because our previous work showed that 
the major impurities are in the form of solid platelets floating on the surface. 
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Thus without the video picture we cannot tell If the 3 ym diameter spot 
analyzed by E«S. Is typical of the entire surface. We have nevertheless 
proceeded vP.h measurements to gain experience and explore the capabilities 
of the device. 

Experl ce ntal Observations 
a. llie Gal' ium Surface 

As we h ve previously reported, (1) the surface of a ne^iy formed 
sessile drop of gallium Is covered by a thin solid oxide layer with a high 
carbon content:. Fig. 4 shows the Auger spectrum of such a drop. In addition 
to oxygen and carbon, this spectrum shows some chlorine which comes from the 
HCl used n the preparation of the drop to dissolve the surface oxide. After 
heating to 410 ^C, the oxygen and chlorine lines are gone, but the carbon is 
undiminish *d as seen in Fig. 5a. A brief sputtering (15 min.) at a beam 
current of 2 x ^0 ^ ampe^ s eliminates the carbon, as seen in Fig. 5b. The 
drop, however, is not clean yet because carbon reappears on the surface after 
a few hours. Several sputtering and equilibration cycles are necessary 
before the carbon is effectively eliminated from the drop. Eventually, no 
carbon Is detected in the spectrum taken at the apex after sitting for many 
days. Although W 2 carrot yet examine for particulates with a video imager, 
we believe the surface to be essentially clean after this procedure. 

Analysis of areas on the drop support this conclusion. 

The r' .uctls>n of carbon concentration in the drop by sequential 
sputtering and equilibration is probably evidence for volume mixing of 
particulates by fl' id flow during Ion bombardment. The Implngent Ions 


remove some carbon from the surface and drive some Into the gallium where it 

Is mixed throughout the drop by convection. Carbon is much less dense than 

gallium so buoyancy win drive these small particles back to the drop apex. 

*4 

A calculation using Stokes equation and a particle diameter of 10 cm gives 

velocities which would clear Che drop volume of carbon in several hours. This 

particle diameter is in reasonable accord with measurements from video images 

2 

of che surface. Since the velocity varies as r > however, it is quite sensitive 
to the diameter chosen, so this calculation only establishes the feasibility 
of this process as an explanation for the evolution of the surface carbon 
concentration. 

There are a number of possible mechanisms which may account for the fluid 
flow during ion bombardment: 1) temperature gradients at the surface, 2) 

electrostatic forces due to surface charging, 3) surface tension gradients due 
to concentration variation of impurities, 4) surface tension gradients due to 
implanted argon concentration gradients, 5) surface density gradients due to 
Implanted argon, 6) momentum transfer from the ion beam. The first of these 
seems imlikely for at least two reasons: a) an ion beam at the low power used 
here does not produce a change in the drop temperature as measured by a 
thermocouple immersed in the drop, b) the impact of the electron beam by 
itself does not produce any mass flow even though its total power is com- 

_3 

parable to that of the ion beam (10 watt). The second explanation also 
seems unlikely since carbon is conductive, as is, of course, the metal 
surface. Surface tension gradients due to Impurity concentration gradients 
should not be present for a surface whiC'. is clean except for a few isolated 
particulates. The fourth and fifth mechanisms are ruled out because they 
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require a significant argon concentration In the gallium and we have never 
been able to detect any argon in the Auger spectra during sputtering. 

The most likely explanation for the fluid flow generated by argon ion 
bombardment seems to be momentum transfer from the ion beam to the liquid. 
Although the magnitude of this momemtum transfer is small, other fluid flow 
phenomena such as thermocapillary convection are also driven by very small 
forces. The rigorous calculation of the flow generated in a sessile drop, or 
any fluid body, by a localized beam of ions would be a major research effort. 
A very approximate calculation, however, yields flow velocities which suggest 

this mechanism to be the one of importance. The ion beam exerts a pressure 

-22 -12 
of about 4 X 10 dvnes/cm over its impact area of 10 cm . If we consider 

the sessile drop surface under the ion beam to be a plane fluid-vapor inter- 
face and neglect fluid flow in tne vapor, the Navier-Stokes equation for an 
incompressible fluid gives 


P 

v 



3V 


2n, 


L dz 


( 1 ) 


for the velocity gradient normal to the surface. Here and are the 

pressures in the beam and in the liquid at the surface, is the velocity 

% -2 

normal to t^ie surface a:.d is the viscosity. For gallium, « 2 x 10 
poise, and using the approximate ion b<*am pressure calculated above gives a 
velocity gradient of 1 s Assuming a linear velocity function normal to 
the surface and a liquid depth of 1 cm yields a normal velocity at the top 
surface of 1 cm/sec. Although this is admittedly a very crude approximation, 
the velocity calculated is large enough to suggest the beam pressure is 
responsible for the fluid flow. 
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Carbon la normally found along with oxygen In the Auger spectra of all 
surfaces before cleaning. Thus, It Is not surprising that we find carbon on 
our gallium surfaces. Gallium drops contained in cups of high density 
graphite show Intolerable carbon concentrations and It has proved impossible 
to eliminate carbon when using them. We presume that the surface of the cup 
is eroding at a rate tdiich keeps the gallium contaminated with carbon 
particles. The vitreous graphite cups used In the Auger work reported last 
year had a much harder surface than the present cups and the gallium could be 
cleaned by repeated sputtering. Even with quartz cups, however, carbon is 
found on the surface of a freshly formed drop as shown in the spectrum of 
Fig. 4. He suspected that the carbon might be originating In the alcohol 
which Is used In the preparation of the gallium before filling the syringe 
for cup loading. However, similar carbon concentrations were found when 
distilled water was used as the carrier fluid. In these experiments, a fresh 
gallium Ingot was used and all tools and glassware were degreased with 
acetone cuid rinsed with distilled water prior to use. 

In a further attempt to Isolate the source of the carbon, a gallium drop 
wh-*u.i had been repeatedly sputtered and showed no indications of carbon after 
e<failx'l itlng for three days was exposed to air by opening a valve. On re- 
evacuation, the Auger spectrum showed a very large oxygen line but no carbon. 
Heating to 500 ”C eliminated the oxygen and the spectrum looked as clean as 
that of a sputtered surface. After equilibrating for two days, however, a 
strong carbon line was measured. It does not seem likely that this carbon 
came from the bulk of the gallium. Rather we suspect that it may have 
originated in the hot filaments of the ionization gauge or the heaters. The 
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residual oxygen pressure due to opening the vacuum system was much higher 

.9 

than normal for several days as the system pumped down to the 10 torr 
range. This oxygen may have reacted with the carbon In the hot filaments to 
generate CO which accumulated on the gallium surface. However, the surface 
showed no oxygen after this procedure so we would have to argue that the CO 
decomposed on the cool gallium surface. This does occur on other metal 
surfaces. The origin of the carbon found on the surfaces of gallium Is still 
unknown. 

b. The Surface Tension of Pure Gallium 

The surface tension of a gallium drop which had been cleaned by sput- 
tering Is shown In Fig. 6. A number of sputtering cycles reduced the carbon 
concentration to the degree that none could be detected twelve hours after 
sputtering. The surface tension values were calculated from photographs of 
the drop profile using the Bashforth and Adans technique. A linear re- 
gression analysis of the data gave the following relationship: 

Y - 724.1 - .0721 T (“C) (2) 

We also show In Fig. 6 the nearly quadratic temperature dependence which we 
found previously for gallium idilch was cleaned by heating In vacuum and in 
hydrogen atmospheres with graphite and quartz cups. This nonlinear re- 
lationship Is In good agreement with three other measurements 
Our present data Is in fair agreement with one of the two linear relation- 
ships found previously for gallium^ and shorn In Fig. 6. Since the Auger 
spectrum shows the surfaces used In the present measurements are clean, we 
believe the nonlinear surface tension temperature dependence found in our 
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Previous work and also In other laboratories may be due to solid carbon 
platelets distorting the sessile drops. These platelets accumulate at the 
apex of the drop, a uniquely important area for surface tension calculation 
because it is the origin from which all the altitudes are measured. Only 
very small distortions would be required to produce the maximum low tempera- 
Cure discrepancy between our two measurements of lOmJ/m at 3U °C, a 
difference of less chan 1.5%. This explanation can only be proved by 
measuring a quadratically varying surface tension for a thermally cleaned 
drop, characterizing the carbon concentration, sputtering to eliminate the 
carbon, then remeasuring the surface tension. We have not done this. 

c. Gallium-Tin 

The first liquid Ga-Sn sample studied in Che Auger spectrometer was a 
gallium drop containing 1.9 atomic percent tin. This alloy was chosen 
because measurements had shown a large reduction in its surface tension below 
that of pure gallium suggesting that tin was strongly adsorbed at the metal- 
vacuum Interface. The sample, however, was contained in a graphite cup and 
we were unable to completely eliminate carbon from the drop by sputtering, as 
discussed earlier. The amplitudes of the gallium and tin Auger electron 
lines are sensitive Co the presence of carbon so that the estimation of 
surface concentration is inaccurate if the surface is not clean. This 
problem was circiimvanted by recording the Auger spectrum after sputtering 
when the surface was free of carbon. Fig. 7 shows the Auger spectrum of the 
alloy drop at 30 and 450 *C. '*he amplitude of the gallium line varies 
little while that of the tin line decreases strongly as temperature Increases. 
Since the amplitude of Che line is roughly proportional to the concentration 
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of the species, this Implies adsorption weakens with increasing temperature. 

This is in accord with observations of surface adsorption in other liquid 
alloys. 

Thu concentration of Sn at the surface of the drop was estimated by 

comparing the amplitude of the dominant tin line at 430 e.v. with that from 

a pure tin sample using the same spectrometer settings. This procedure is 

not straightforward, however, because the detected Auger electrons originate 

from significant depths and -.lot solely from the surface layer. The primary 

electron beam (200 e.v. for our measurements) penetrates deeply into the 

sample. The Auger electrons generated by the beam are exponentially 

attenuated with distance from their origins. Thus when Auger line amplitudes 

from surface layers are compared with those from bulk standards, this 

[81 

attenuation must be taken into account^ . For an incoherent scattering 

model of a uniformly adsorbed layer of Sn, the intensity of the 430 e.v. 

Sn 

Auger electrons, , is given by 


I 


Sn 



1 + (430 e.v.) 

1 + r_ (430 e.v.) 
Sn 



(3) 


Here is the intensity of a pure bulk tin sample, m is the layer thickness 

and L is the attenuation length for 430 e.v. electrons in tin. The r terms 
take into account the Auger electrons which arise from the back&cattered 
primary beam; they vary with energy and atomic species. In the filled circles 

in Fig. 8 we show the surface coverage of tin as a function of temperature 

[91 

calculated with equation 3 using an attenuation length of 4 monolayers • 

The r values for Sn and Ga at 430 e.v. were estiniated to be 0.62 and 0.46 


19 


Is approximately 0.9. 


respectively . Thus the term -rTT^ 

Sn 

We have also measured the surface tension for the gallium-1. 9% Sn drop. 

The data shown in Fig. 9 uHled circles on line 1) were taken after sputtering 
and before carbon redeveloped at the surface. Although carbon may have been 
located on the surface outside of the area analyzed by the electron beam, t'ue 
error so introduced into i.he surface tension values should be small compared 
to the large effects of the tin adsorbed at the surface. We also show in thlu 
figure the surface tension relationship (line 4) we find fo-. pure gr Ilium 
surfaces cleaned by sputtering. For comparison some surface tension values 
for gallium and a galllum-2% tin alloy measured in another laboratory (ref. 6) 
are also presented. The results for the alloys are in fairly good agreement 
except below 100 ®C. 

The depression of the surface tension of the alloy below that of pure 

gallium is produced by the adsorption of tin at the liquid-vacuum Interface. 

This adsorption is strongly temperature dependent and produces a temperature 

dependent surface concentration of tin, as indicated In Fig. 8. This surface 

(Ga) 

concentration of tin can be estimated from the relative adsorption F ' 

on 

For a dilute solution this can be written as 


r (Ga) . _ ^Sn 

Sn RT dX- ^ ' 

Sn 

where Xg^ is the bulk concentration. The coverage in monolayers is found by 
/ \ 

multiplying by the area of a tin molecule. In Fig. 8 we show (open 

circle) the surface coverage of Sn for the gallium alloy drop calculated by 
taking 




alloy 


( 5 ) 


20 




These results can be compared with the surface coverages estimated from the 
Auger line intensities (filled circles in Fig. 8). The surface coverages are 
in good agreement, particularly a. the lower temperatures. Given the uncertainties 
in attenuation lengths and backscattering factors, the agreement may be 
fortuitous. No effort has been made to improve the agreement of the two sets 
of data by adjusting these parameters. The surprising aspect of this data is 
tliat it suggests coverages in excess of a monolayer below 125 *^0. It was 
anticipated that the coverage might reach a monolayer but would then saturate 
and remain constant. Large adjustments to the attenuation length, back- 
scattering factors, and the variation of surface tension with concentration 
would be necessary to reduce the coverage to below a monolayer. 

In conclusion, our initial studies have shown that Auger spectroscopy can 
be successfully applied to liquid metals and alloys. Ionic bombardment 
sputtering seems to be a promising way of producing clean liquid surfaces 
although we have not tested the technique yet for impurities which are surface 
active but remain in solution. The study of segregation at liquid alloy 
surfaces by Auger spectroscopy combined with surface tensxoii measurements is a 
promising area for future work even though attenuation lengths and back- 
scattering factors may not be sufficiently well known for f.ome materials. 
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Photograph (negative) of sessile drop in Auger spectrometer. 



Fig. 4. Auger spectrum of uncleaned Ga drop. 
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Fig. 5. Auger spectrum of Ga drop (a) after heating to A 50 ®C (b) after sputtering. 






The Auger spectrum of Ga-Sn 




M (monolayers) 



I I I I I 
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Fig. 8. The surface coverage in monolayers of Sn as a function of temperature 
as measured from line amplitudes (filled circles) and surface 
tension chai\ges (open circles). 
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Fig* 9* The surface tension of Ga-Sn 1*9 at*X as a function of temperature* 
The filled circles aloig with curve (1) give measurements from the 
current work* The X's which form curve (2) are data reported in 
reference 6* Curves (3) and (4) are results for pure gallium, the 
X*s which form curve (3) being data given in reference 6 and the 
straight line (4) being results from the present work. 
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Task 2 


Convection During Unidirectional Solidification 

S. R* Coriell and R. J. Schaefer 
Metallurgy Division 
Center for Materials Science 

Summary 

The succlnonitr 11 e-ethanol system has been selected for the experimental 

study of solute- induced convective and morphological instabilities • The phase 

diagram of this system has been determined and the measured distribution 

coefficient of 0.044 vas used to perform stability calculations. Several 

-3 -2 

modes of instability are predicted In the composition range of 10 to 10 
weight percent ethanol and the solidification velocity range of 0.5 to 5 pm/sec. 
Experimental measurements of convective flow have used the tracking of small 
neutrally buoyant particles, and have to this point been devoted to measuring 
the background of thermally induced convection in pure succlnonitrile. Flow 
rates of 1 to 5 pm/ sec were produced in the present apparatus from thermal 
convection alone without solidification whereas convective flows induced by 
solidification processes were 15 pm/sec or more. 


Preceding page blank 
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Introduction 


During the solidification of a material containing more than one 
component, density differences due to chemical inhomogeneities of the 
liquid phase are among the several driving forces for convective flow, 
in the presence of a gravitational field* The onset of coupled con- 
vective and constitutional interfacial instabilities during the directional 
solidification of a single phase binary alloy at constant velocity 
vertically upwards was previously treated by a linear stability analysis 
and specific calculations were made for physical properties appropriate 
to the solidification of lead containing tin [1]« For experimental 
verification of this type of calculation, and to reveal the nature of 
the convective flow once it is Initiated, it is necessary to study the 
solidification of transparent materials. Such materials have physical 
properties which differ greatly from those of metals, but similar 
instabilities should occur. We report here calculations and experiments 
on the solidification of succinonitrile containing ethanol. For the 
lead-tin system^ we give some additional results in which the convective 
and interfacial instabilities are decoupled and compare these with the 
coupled results. 


Theory 

We recall that the dependence of the perturbed temperature, con- 
centration, fluid velocity, and interface shape on the horizontal 

coordinates x and y and the time t is of the form exp [at+iCwx+wy)], 

X y 

where a> and o) are spatial frequencies. In general, the time constant 
X y 

o is complex, i.e., a ■ a + io.. If o >0 for any values of caj and u) , 
-r f f r i r ^ X y 
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the system is unstable. If a <0 for all values of u> and u) , the 
^ X my 

system is stable. We refer to Instabilities corresponding to large 

2 2 % 

values of u) « (u> + u> ) (small wavelengths, X » Ztt/co) as morphological 

^ y 

(or interfacial) instabilities and to instabilities corresponding to 
small values of o) (large wavelengths) as convective instabilities. 

The properties of the succlnonitrile-ethanol system are given in 
Table I. The solute diffusion coefficient, D, in the liquid and the 
variation of liquid density, p, with solute concentration, c, have not 
been measured. Estimated values of these have been used in the calcu- 
lations, viz., D * 1.0 (10 ^) cm^/s and e - (3p/3c)/p * 3.07 (10 ^) 


.-1 


For a temperature gradient G in the liquid of 10 K/cm, the 

ii 


(wt.%) 

critical concentrations of ethanol in succinonitrile above which in- 
stability occurs are shown in Fig. 1. The results are similar to the 

previous calculations on the lead-tin system, and instability occurs at 

-3 

concentrations of the order of 10 wt.%. For small velocities the 

convective instability occurs at lower concentrations than the mor- 
phological instability. The concentration for convective instability 
increases with velocity while the concentration for morphological 
instability decreases with velocity. The two curves uoss at a velocity 
of about 3.5 pm/s. In the vicinity of this crossover, i.e. , where 
convective and morphological occur at about the same critical concen- 
trations, instability occurs by an oscillatory mode for which Q 

when =» 0. The onset of instability by an oscillatory mode is the 
most striking dlffer'^nce between the results for lead- tin and succinonitrile- 
ethanol. Fig. 1 also shows the dimensionless ratio of the perturbation 
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wavelength A at the om^et of instability to the diffusion boundary layer 
thickness D/V* As previously noted, the wavelength corresponding to 
convective instability is considerably larger than the wavelength corre- 
sponding to interfacial instability. At small velocities, the convective 
wavelengths are quite large, e.g. , A = 5.7 cm for V « 0.5 um/s. 

The effect of the temperature gradient in the 1 quid on the 

critical concentration for instability is shown in Fig. 2. As long as 
the mode of instability does not change, the critical concentration 
increases with Increasing temperature gradient as would be expected 
since the temperature gradient is a stabilizing influence. However, as 
illustrated by the V 3.0 ym/s data in Fig. 2, the critical concentration 
at G^ * 10 K/cm is larger than the critical concentration at 20 K/cm; 
the mode of instability changes from oscillatory convective (a^ ^ 0) to 
convective (o^ * 0) as the gradient increases. For V = 4.0 ym/s, the 
IT de of instability changes from morphological to oscillatory convective 
to convective as the temperature gradient increases from 10 to 40 K/cm. 

These changes in the mode of instability are illustrated in detail 

in Figs. 3-6, which give the concentration at which = 0 as a function 

of the spatial frequency o) of the sinusoidal perturbation for * 10 

*v/cm and V - 2,5, 3.0, 3.5, and 4.0 ym/s, respectively. The solid lines 

correspond to * 0, while the dashed lines correspond to ^ 0; the 

values of are indicated by the lower curves in the figures. At V = 

2.5 ym/s, the smallest concentration corresponds to a non-oscillatory 

convective mode whereas at V * 3.0 and 3.5 ym/s an oscillatory convective 

-2 -3 -1 

mode occurs with in the range 10 to 10 s . At V * 4.0 ym/s, 

-3 

instability first occurs by a morphological mode with c^ * 4.66 (10 ) wt.% 
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V 


and u) = 140 cm ^ (not shown in the figure). 

Since the liquid diffusion coefficient D has not been measured, some 

calculations were carried out with a different value of D. For D = 5 
—6 2 

(10~ ) cm /s and G = 10.0 K/cm, tb« critical concentrations are 2.09 

Xj 

I 5.06 (10“^), and 3.46 (10“^) wt.% for V * 0.5, 1.0, and 2 0 ym/s, 

I 

* respectively. The highest velocity corresponds to an oscillatory convective 

i 

instability with non'-osclllatory convective instability occuri.xng at the 

-5 2 

lower velocities. For comparison, for D = 1.0 (10 ) cm /s and = 

-4 ~4 

10.0 K/cm, the critical concentrations are 2.63 (10 ), 4.58 (10 ), 

_3 

and 1.20 (10 ) wt.% for V = 0.5, 1.0, and 2.0 ym/s, respectively. 

Thus, the critical concentration curves (as function of velocity) for 

the two different diffusion coefficients cross each other. The critical 

concentration for morphological instability is proportional to D for 

small velocities, and hence decreases with decreasing D. 

In order to obtain an estimate of the rate at which instabilities 

develop, linear stability calculations of o for a fixed concentration of 

-5 2 

ethancl were carried out. For D - 1 (10 ) cm /s, G = 10 K/cm, and 

-2 

c^ = 10 wt.%. and crystal growth vclocites of 1.0, 2.0, and 4.0 um/s, 

-2 

the maximum values of o for convective instability are 4.9 (10 ), 

-2 -3 -1 

3.7 (10 ), and 6.6 (10 ) s at spatial frequencies of 22, 24, and 

28 cm respectively. For V = 2.0 and 4,0 ym/s, the maximum values of 

-2 ^1 

0^ for morphological instability are 1.1 (10 ) and 0.12 s at spatial 

frequencies of 440 and 950 cm respectively. For V = 1.0 ym/s, the 
system is morphologically stable, i.e. , the maximum value of for 

morphological instability is negative. For example, at V = 1.0 ym/s, 
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the amplitude of the convective instabilit> will increase by a factor of 
ten in 50 s; during this time the crystal will grow 0.005 cm. 

In order to gain some insight into the interaction between convective 
and morphological instabxlities, we have modified certain parameters ^n 
the calculations so that only one type of instabilxty can occur. Physical 
parameters appropriate to the solidification of lead containing tin were 
used with V = 40.0 um/s and = 200 K/cm, If the gravitational acceleration 
g is identically zero, then convection does not occur and only morphological 
instabilities are possible. The critical concentration of tin as a function 
of the spatial frequency u of a sinusoidal perturbation with g = 0 is 
shown in Fig. 7; the minimum value of c^ is 0.215 wt.% at w = 690 cm By 
requiring that the solid-liquid interface not deform, we can eliminate 
the morphological instability and allow only convective instabilities. 

The boundary conditions at the solid-liquid interface are based on con- 
servation of enthalpy and solute and a relationship between freezing 
point and solute concentration. In order to simultaneously satisfy 
these three boundary condifLons, the solid-liquid interface must be free 
to deform. If we require a planar interface, one of the boundary con- 
ditions can not be satisfied. We have required that the interface be 
rigid (planar) and satisfied the coiiservation laws but ignored the 
relationship between freezing point and solute concentration. The critical 
concentration as a function of w for this rigid interface model is shown 
in Fig. 8; the minimum value of c is 0.15 wt.% at w = 45 cm The 

oo 

2 

results for the coupled problem, i.e. , g = 980 cm/s and Interface free 

to deform is shown in Fig. 9. The minimum value of c is 0.136 wt.% at 

" 00 
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0 ) » 45 cm ^ for convective instability while the minimum for the morphological 
Instability is essentially identical to that computed with g » 0. 

While there is a 25Z difference in the critical concentrations for 
convective instability between the rigid and free interface, the most 
strileing difference between Fig. 9 and Figs. 7 and 8 are the appearance 
of oscillatory modes (a^ ^ 0) in the coupled problem. These oscillatory 
mooes are apparently the result of the coupling between the convective 
and morphological instabilities, and would not be expected in the absence 
of this coupling. Oscillatory modes do arise in double diffusive con- 
vection, but in the simple models which can be treated analytically, 
only when the temperature field is destabilizing and the solute field 
stabilizing, which is not the case treated here. The calculations for 
the succinonitrile-ethanol system have indicated that the onset of 
instability can occur by an oscillatory mode. 

The main theor^*"ical effort of this task is now directed toward 
developing mathematical models cf che fluid flow and resultant solute 
segregation under conditions for which the linear analysis predicts 
convective instability; this research is being carried out in collaboration 
with R. G. Rehm of the Mathematical Analysis Division. Numerical 
algorithms are being developed and implemented to solve the time dependent 
partial differentiax equations for fluid flow and heat and mass transfer 
in two spatial dimensions with a planar solid-liquid interface. 

During the past year, a paper "Effect of Gravity on Coupled Convective 
and Interfacial Instabilities During Directional Solidification," by 
S. R. Coriell, M. R. Cordes, W. J. Boettinger and R. F. Sekerka, was 
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presentisd at the XXIIl COSPAR seetlng and vlll be publiahcKl in the 
proceedings. An article, ”Solutal Convection Induced Hacrosegregation 
and the Dendrite to Composite Transition In Off-Eutectic Alloys/' by W. 
J. Boettinger, F. S. Biancaniello, and S. R. Coriell, which describes 
previous work of this project, has been published in Metallurgical 
Transactions [ 2 ] . 

Experimental Procedures 

Because no methods are available for the detailed measurement of 
fluid flow patterns in liquid metals, a transparent system is used. 
Succinonitrile is chosen as the major constituent because of its metal- 
like solidification behavior and well-characterized physical properties 
(Table I). Extremely high purity can be attained in this substance by 
multiple distillation and zone refining [3J. Unfortunately, it has a 
rather low density, so the choice of solutes having lower density is 
extremely limited. Ethanol has been selected for this purpose on the 
basis of its low density, relatively high boiling point compared to 
other low density materials, and similarity of molecular structure to 
succinonitrile * 

For input to the stability calculations we need to know the 
succinonitrile-rich end of the succinonitrile-ethanol phase diagram. 

This was determined by observations of samples immersed in a stirred 
water thermostatic bath* Temperatures were measured with precision 
thermometers on which temperature differences as small as 0*002 K can be 
read* The temperatures are not known with comparable precision on an 
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absolute scale, but only Che accurate teaperature differences are liq>ortant 
for this experlaent. The 8ucclnonltrlle->eChanol systea has a liquid- 
phase alsclblllCy gap, as does succlnonltrlle with several other alcohols. 
The llquldus curve (Fig. 10) was deteralned froa saaples weighed out and 
sealed -aider air. A aonotectlc point was found at approxiaately 11 *C 
and 23 wt.Z ethanol. The odscibility gap is rather low but was not 
studied In detail because the region of interest lies at very low con- 
centrations of ethanol. The llquldus curve is essentially linear at 
concentrations of ethanol of less than 6 wt.Z, and its slope in this 
region is 3.6 K/vt.Z. 

The solidus curve is very steep and its detemlnatlon could only be 
carried out with saaples held under vacuum, because dissolved air acts 
as a significant impurity. Such samples are difficult to prepare with 
accurate compositions because the vapor pressure of ethanol is high at 
the melting point of succlnonltrlle. Therefore, the apparatus shown in 
Fig. 11 is employed. Zone-refined succlnonltrlle Is first loaded Into 
the larger chamber, melted and frozen several times with active pumping 
to remove all dissolved air, and this side of the chamber Is then 
hermetically sealed. The mass of the succlnonltrlle Is determined by 
weighing. Ethanol Is then loaded Into the capillary chamber and part of 
It Is frozen unidlrectlonally by lowering the capillary into a bath of 
liquid nitrogen. This side of the chamber Is then evacuated and sealed. 

Additions of ethanol to the succlnonltrlle are carried out by 
opening the teflon valve and gently applying heat to the top of the 
ethanol column to distill over a small amount. An estimate of the 
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amount transferred is obtained by noting the drop in the ethanol level in the 
capillary. A more accurate measure of the concentration of ethanol is obtained 
from the measured liquldus temperature (the tanperature at which only a minute 
fraction of the sample remains solid, and the solid is neither melting nor 
freezing). The concentrations determined by these two methods agreed well. 

The solidus temperatures were determined as the lowest temperature at which 
melting was first observed (at grain boundaries) after the sample had been 
honogenized by rapid freezing. The results of these experiments are shown in 
Fig. 12, and the slope of the solidus curve is found to be 81 K/wt.%. 

The distribution coefficient is a key parameter in the theory of solutal 
convection since this coefficient determines the proportion of solute rejected 
at a solidifying Interface. With the above data, this parameter can be cal- 
culated from the ratio of the liquldus slope to the solidus slope at small 

solute concentrations to be k » 0.044. 

o 

To investigate the convective and interfacial instability of succinonitrile- 
ethanol solutions, the calculations have indicated that one should grow crystals 
in the solidification velocity range 0.5 to 10 pm/s. For these experiments 
the sample is sealed in a round pyrex tube, 19 mm in diameter. This tube is 
dram downward at the desired velocity through the apparatus shown in Fig. 13. 
This apparatus contains an upper water jacket, maintained at a temperature 
above the melting point of succinonitrile, and a lower Jacket maintained below 
the melting point. Between these reservoirs, the round sample tube is enclosed 
in a fluid-filled Jacket with flat glass walls. The fluid in this Jacket is a 
mixture of water and ethylene glycol having a refractive index such that the 
optical distortion of the material within the round sample tube is minimized 
[3]. The upper and lower water Jackets are maintained by circulating controllers 
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at temperatures constant within + 0.001 K. The apparatus is heavily Insulated 
and the sample Is viewed In the region of the solld^llquld Interface through 
small double-walled glass windows. 

Convection is measured by photographic tracking of the motion of small 
particles suspended In the liquid. Laser Doppler Veloclmetry (LDV) was con- 
sidered as a measurement technique but rejected because at the extremely slow 
flow velocities (down to 1 pm/s) expected in these experiments, the LDV technique 
would be overly sensitive to mechanical vibration and electronic noise, and 
because it is more important for these experiments to determine the overall 
pattern of the flow than to obtain statistical information on velocities at a 
fixed point in space, as one normally does with the LDV method. 

The particles used for tracking the fluid flow are polyvinyltoluene 
spheres 2.02 urn in diameter. The small difference between the density of 
these spheres and that of the liquid succinonitrile, together with their small 
size, results in a rate of gravitational precipitation calculated from Stokes 
Law to be 0.05 iim/s, so that the spheres follow almost perfectly the flow of 
the liquid. The particles are readily visible through a microscope, when 
illuminated by low-angle dark field illumination. Photographs are recorded 
with exposure times usually between 30 seconds and 2 minutes, to give adequate 
travel distance for flow rate measurements. The particles are illuminated by 
a microscope light, with a heat-absorbing solution in the optical path to 
prevent radiant heating of the specimen. 

Experimental Results 

During slow solidification, the particles are pushed by the solid/ 
liquid interface, often accumulating in the small grooves which are 
present where grain boundaries intersect the solid/liquid interface. 
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Vfhen the partitles accumulate Into large clusters, or when the interface 

velocity becomes rapid, the particles are Incorporated into the solid. 

No quantitative measurements of this effect have been recorded, and 

it does not interfere with the convection measuremencs. Convection 

measurenients to date have used a sample of ultra high purity sucl I nonltrlle 

in which all convection must be attributed to thermal effects. This 

sample was prepared by the same techniques which were previously found 

[3] to reduce the level of impurities to approximately 5 x 10 ^ to 5 x 
—8 

10 mole fractions The melting and solidification of this sample shows 
no indication of ljrpurities» such as grain boundary melting or mor- 
phological instability. The objective of the measurements with the high 
purity sample is to establish a baseline of residual thermally- induced 
convection to which convection in ethanol-doped samples can be compared. 

When the temperature reservoirs of the crystal growth apparatus -re 
adjusted to give a gradient of about 10 K/cm in the liquid, and the 
sample tube is held stationary (no crystal growth), the solid- liquid 
interface is planar and horizontal, and convective flow rates in the 
liquid are generally foimd to be in the range of 1-5 \m/^. The flow 
pattern is somewhat erratic. An occasional source of higher flow rates 
under these conditions was seen when a large void, which had originated 
as a shrinkage cavity when the sample was originally solidified, migrated 
vertically upward through the solid under the influence of the temperature 
gradient. Shortly before this void reached the solid-liquid interface 
(where it collapsed), it was seen to have an influence on the convective 
flow field such that there was a visible upward flow of the liquid above 


44 ‘ 



the void* The thermal distortion which produced this flov» however » was 
not sufficient to produce a measurable distortion of the solld*liquid 
interface. 

When a stationary sample is suddenly repositioned downward, crystal 
growth starts almost iinnedlately, with the most rapid growth along the 
tube walls. As expected, convective flow is greatly accelerated. 

Several hours are required before the interface returns to its original 
position and the convective flow again diminishes to velocities char- 
acteristic of stationary samples. When the sample was displaced 1 cm, 
the flow rate along the sample axis was 20 um/s at 40 minutes after the 
displacement. The flow pattern was a steady flow upward along the 
sample axis and downward along the walls. 

When the mechanism for pulling the sample downward at constant 
velocity is activated, the solid-liquid interface moves downward and 
eventually reaches a new steady-state position. Typically 2-3 hours are 
required before the interface position becomes effectively stationary. 

In addition, the interface changes from the planar horizontal configuration 
to an upwardly concave shape. For a sample which was drawn downward at 
a velocity of 2 Mm/s, convective flow rates were found to increase as 
the interface became more concave, so that one hour after the tracking 
mechanism was started, flow rates were 15-20 ym/s. 

Although the convective flow during constant velocity growth is too 
small to be readily visible to the eye when viewed through the micro- 
scope, it is desirable to modify the apparatus such that this purely 
thermal effect is minimized. To this end, experiments will be carried 
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out with Auxiliary heatera in the index-matching flu:!'! ao that the 
diacontinuity of temperature gradient induced by the latent heat of 
fualon can be approximated in this fluid. Experiments with ethanol- 
doped samples will be connenced to determine the relative magnitude of 
the flow rates in pure and doped samples. Simultaneously improvements 
in sample illumination will be sought which can reduce the intensity of 
light scattered from extraneous sources such as the chamber walls and 
thus allow a much wider range of velocity measurements from a single 
exposure. 
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Table I. Properties of Succlnonltrile 


Property 

Symbol 

Value 

Ref. 

Molecular Weight 

w 

80.092 

4 

Density of Solid 

^s 

1016kg/m^** *** **** 

4* 

Density of Liquid 


988kg/m^** 

i 

5 

Thermal Expansivity of Su^ld 

U 

s 

-5.6x10 /K 

6 

Thermal Expansivity of Liquid 


-8.1xlO”^/K 

5* 

Shear Viscosity 

n 

2.6xl0‘^Pa-s** 

5 

Kinematic Viscosity 

V 

2.6 xlO ^m^/s** 


Surface Tension (Liquid-Vapor) 

^£V 

46. 78mJ/m^ 

5 

Surface Tension (Solid-Liquid) 


8. 9mJ/m^ 

3 

Refractive Index (Solid) 

n 

s 

1.4340** 

7 

Refractive Index (Liquid) 


1.4150** 

7 

Equilibrium Temperature (Triple 


331. 23K 

8 

Point) 

Latent Heat of Fusion 

L 

4.70x:0^J/kg 

4 

Entropy of Fusion 

AS 

1. 42x10* j/kgK 

4 

Heat Capacity oi Solid 

Ss 

1913j/kgK** 

4 

Heat Capacity of Liquid 


2000J/kgK** 

4 

Thermal Conductivity of Solid 

k 

s 

0.225j/mKs** 

3 

Thermal Conductivity of Liquid 

s 

0.223j/mKs** 

3 

Thermal Di^fusivity of Solid 

1.16xl0“^m^/s** 


Thermal Diffasivity of Liquid 


1.12xl0'^m^/s** 


Unit Supercooling 

L/Cp^ 

23.5 K 


Distribution Coefficient for 
Ethanol in Succlnonltrile 

k 

o 

0,044 


Liquidus slope for Ethanol 

m n 

3.6K/wt.% 

iciCTkis 

in Succlnonltrile 

Solidus slope for Ethanol in 
in Succinonitrile 

m 

s 

81K/wt.% 



* Indicates value derived from properties given in cited reference. 

Indicates value of property at melting point, 331. 23K. 

*** Indicates property der-*ved from other properties in the table. 

****Indicates properties measured in current work. 
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1.0 


10.0 


V /im/s 

Fig. 1. The critical concentration c^ of ethanol in succinonitrile above 
which instability occurs as a function of the velocity V of 
directional solidification for a temperature gradient in the 
liquid of 10 K/cm. The ratio of the instability wavelength, X, 
to the diffusion boundary layer thickness, D/V, is also shown. 

The solid and dashed lines correspond to interfacial and cca- 
vective instabilities, respectively, while the dotted line denotes 
an oscillatory convective instability. 
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SUCCINONITRLE 

ETHANOL 

V*2.5/xm/s 

6,*10K/cm 


Fig. 3. The concentration of. ethanol in succinonitrile at the onset of 

instability during directional solidification at V = 2.5 ijm/s as 
a function of the spatial frequency to or a slmisoidal perturbation 
The solid curves ma’-fc the onset of non-oscillatory instabilities 
(o *0); whereas the dashed curves mark the onset of oscillatory 
instabilities (the value of o. is given on the right hand side). 





10^ Coowi% 



Cl/cm 


Fig. 4. The concentration of ethanol in succinonitrile at the onset of 

instability during directional solidification at V = 3.0 ym/s as 
a function of the spatial trequency w of a sinusoidal perturbation. 
The solid curves mark the onset of non-oscillatory instabilities 
(o.=0); whereas the dashed curves mark the onset of oscillatory 
instabilities (the value of is given on the right hand side). 
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SUCaNONITRILE-ETHANOL 
V*3.5 ^m/8 

6i *K> K/cm 


5 o 


The concentration of ethanol in succinunitrile at the onset of 
instability during directional solidification at V = 3.5 pm/s as 
a function of the spatial frequency to of a sinusoidal perturbation. 
The solid curves mark the onset of non-oscillatory instabilities 
(o^=0) ; whereas the dashed curves ‘mark the onset of oscillatory 
instabilities (the value of a. is given on the right hand side). 






IO®CmWt% 



CJcrn'* 


Fig. 6. The concentiation of ethanol in succinonitrile at the onset of 

instability during directional solidification at V = 4.0 ym/s as 
a function of the spatial frequency cu of a sinusoidal perturbation. 
The solid curves mark the on^^et of non- os dilatory instabilities 
(a^=0); whereas the dashed curves mark the onset of oscillatory 
instabilities (the value of o. is given on the rxght hand s^de). 


40 ft.m/% 
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The concentration of tin in lead at the onset of instability during directional 
solidification at V = 40.0 pm/s as a function of the spatial frequency u of a 
sinusoidal perturbation. The gravitational acceleration g « 0 so that only 
morphological instability can occur. 




RIGID INTERFACE 
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Fig. 8. The concentration of tin in lead at the onset of instability during directional 
solidification at V = AO.O vim/s as a function of the spatial frequency . of a 
sinusoidal perturbation. The sol id- liquid interface is rigid (planar) so lliat 
• )iily lonvective instability can occur. 





Teflon Valve 


r\ 



Succinonitrlle Ethanol 

chamber Capillary 


Fig. 11. Apparatus for the preparation of succinonitrile-ethanol mixtures 
under vacuum. 
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WBGHT PERCENT ETHANOL 


Fig. 12. Detail of succlnonitrile-ethanol phase diagram at very small 

concentrations of ethanol, showing data points used to determine 
the solidus slope. 
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fig. 13. Unidirectional crystal growth apparatus. The sample is drawn 

downwards to produce crystal growtJi. Thermal insulation surrounding 
the water jackets and sample tubes is not shown. 
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Task 3 

Measurement of High Temperature Thermophysical Properties 
of Tungsten Liquid and Solid 

D. W. Bonnell 

Chemical Stability and Corrosion Division 
Center for Materials Science 

SUMMARY 

The combined General Electric Advanced Applications Laboratory 
(GE) , Rice University (RICE), National Bureau of Standards (NBS) effort 
has reached the end of its second phase. This phase was intended to 
complete the physical connection of the RICE calorimetric system to the 
GE levitation facility, and establish the limits of the system for the 
acquisition of solid and liquid enthalpy increment data for tungsten in 
an actual measurement effort. These limits have been established, 
indicating that data are obtainable on solid tungsten up to the melting 
point with a precision comparable to prior high temperature levitation 
studies. It appears, however, that while liquid data are obtainable, 
the reliability of drop successes is so low that it will be necessary to 
completely automate the calorimeter system in order to accommodate the 
extremely unreliable holding characteristics of the levitation system. 
This unreliability is exacerbated by the nature of the calorimetric 
measurement requirement . 

A number of necessary improvements in the general system have been 
made, particularly in the temperature measurement area, which allow 
calorimetric quality heat content measurements to be related to an 
accurate temperature. These improvements are expected to form a model 
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which GE will use in a variety of areas, and will hopefully provide the 
basis for a generally available very high tenperature imaging pyronetric 
device. 

Finally, this second phase effort helped establish in detail the 
precise requirements and degree of control necessary (1) to provide a 
major improvement in the ability of this experimental effort to yield 
accurate data and (2) tc act as a model for future efforts in high 
temperature calorimetry, such as space based experiments where only a 
single operator will be available. Such an operator will, of necessity, 
need to depend on the instrumentation to control details and provide 
appropriate feedback, in order to be free to pursue the goals of the 
experiment, without being lost in the simple achievement of operation of 
the experiment. 



INTRODUCTION 


The primary aim of this task has been to evaluate experimental 
procedures used in the interaction between the General Electric Advanced 
Applications Laboratory (GE) and the Rice University (RICE) to measure 
the high temperature enthalpy increments of liquid and solid tungsten. 

GE has demonstrated a unique apparatus capable of levitating molten 
tungsten (1] and RICE has wide experience in high temperature enthalpy 
increment measurements, being particularly successful in the use of 
drop-type isoperibol calorimeters with levitation heating (2;3a,b;4]. 

It seemed a natural cooperative effort for RICE to provide the calorime- 
tric apparatus and GE to provide the levitation facility in order to 
measure the high temperature enthalpy function of tungsten. 

The results of this research are of great theoretical and engineering 
interest. From the scientific viewpoint, the unique location of tungsten 
at the upper extreme of the metal and element melting point scale should 
provide a key part in any extrapolation/interpolation procedure. Models, 
such as th( well-known Tamman (5] rule and periodic table correlations [6], 
are currently used for extrapolation. A definitive value for the heat of 
fusion of this element will provide a valuable test of such uiodels, perhaps 
indicating a significant deviation in such predictions based on lower 
melting refractory metals. Questions concerning the liquid state heat 
capacity function, i.e., whether it is constant as would be suggested 
from studies of lower melting refractory metals (see, for example, (7]) 
and whether current estiro-tes (8;9] of the discontinuity from Cp(s) to 
Cp(l) are even reasonable, are crucial in testing theoretical models 
and may provide clues to establishing new parameters in existing 



or new models. The engineering applications beg^u with the tact that, 
without knowledge of the heat of fusion and liquid enthalpy function, 
melting, forming, and casting equipment must of necessity be overdesigned 
to allow a satisfactory safety/performance margin. These engineering 
effects pervade all high temperature uses of tungsten*containing materials. 

In the course of this work, it has een necessary to be present 
during much of the actual experimental work at GE, and both SICE and GE 
consult this author concerning the resolution of problems at all phases 
of the experimental process, from initial design to final data reduction. 
All modifications to the RICE calorimetric system were cesig*.d by <his 
author, and a large amount of available research time has had to be 
devoted .o the adaptation of the GE observation system to be an accurate 
pyrometiic device. 

Liquid holding capability has proved to be the current unresolved 
problem. The solution seems to require a totally automated calorimeter 
which can be set to detect a falling sample, actuate its gates in real 
time, and acquire the sample, while performing in real time, the prelim- 
inary calculations of data reduction necessary to establish the ready 
state of the calorimeter, M.e true heat influx, and when the calorimeter 
again reaches steady state. Such automation, when coupled to an automatic 
version of the now Ui/er;.tional GE pyrometric system, will provide the 
basis for future experiment, perhaps in space, where real time data 
acquisition is not a luxury, but a necessity. 



EXPERIMENTAL PROCEDURE 


The current experimental configuration is shown schematically in 
Figure 1 (a,b). The sample is supported prior to levitation ::n a rovable 
pedestal, part of a lava and tantalum platform on rails which can be 
laterally moved clear of t^e coil bottom during levitation. This platform 
can be used to load successive spheroidal samples, but has only marginal 
ability to select the order of introduction among a group of samples. 

The tungsten radiation gat*» below the pedestal assembly must be manually 
rotated clear of the drop path to obtain a successful drop. A micro* 
switch on the radiation gate actuator detects the opening of the gate, 
turns off the electron beam, and turns down the RF power to the levita- 
tion coil to drop the sample under control. 

The work coil has been recently been redesigned. The original 
design was a lower-coil-only conical design where the molybde.ium "shade’* 
ring actea as a dock for field closure [10;11). On our recommendation, 
because of severe instability and the inability to melt, the current 
coil design depicted in Figure la. was implemented by GE. This coil 
style is a conical modification of the coil type used for liquid platinium 
cort^inment at RICE [4], incorporating an extra bottOBC turn to reduce 
coil leakage. The requirement to pass the physically targe (~ 1 cm) 
spherical samples and the large pedestal results in a bottom coil throat 
inside diameter of ~ 1.3 cm. This large bott . opening still exaggerates 
the pcr:sibility of sample loss due to drip»j,ii^, as was the case for the 
original coil, and it is certain that more work needs to be done by GE 
in this area. The incorporation of an upper coil turn does, however, 
ailov’ higher coil power to be used to improve stability and lessen the 



probability of coil dripping. The shade ring over the coil has been 
retained to continue to serve as electron beam s^^ll protection for the 
work coil. 

The levitation coil is driven by a 25 kW, 450 kHz General Electric 
induction furnace. The vacuum feedthrough to che coil is a GE coaxial 
design; the entire transmission line including inductive coupling (matching 
transformer) allows power factors > 0.8. 

The levitation coil is not designed to provide sufficient heating to 
melt the sample. The maximum temperature reached by induction power 
alone is approximately ?&Q0 K. A simple calculation for the 1 cm diameter 
samples of this work using the Stefan-Boltzmau Law, i.c., 

dH/dt = E^aA(T^) 

(where dH/dt is the rate of radiant energy loss for a body of area A, 
total hemispheric emissivity E.p at a differential temperature A(T) to 
its surroundings, a is the Stefan-Boltzmann constant of proportionality, 
5.67032 X 10 W*cm ^ [12]) gives the 2800 K radiation loss as 

about 310 watts when a value for E^ of 0.28 [13] is chosen for tungsten. 

At ♦'he melting point of tungsten, 3695 K [14], the loss is expected 
to be about 1000 watts. This additional power is provided by electron 
beam (e-beam) heating. In the levitation environment, e-beam heating is 
possible only at temperatures where secondary electron emission is high 
enough to maintain approximate charge neutrality. Tungsten ^s secondary 
emission characteristic ro^ets this criterion even for electron power 
fluxes much in excess of 1 kW [15]. The electron beam gun shown in 
Figure la directs a 35 kV, 150 mA focused electron beam onto the top of 
a levitated specimen. This power is clearly sufficient to achieve 
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■elting even with considerable loss due to charging. A significant 
problen in this form of heating is the localized influx of energy due to 
spot focusing. While the sample is solid, the top to bottom teaiperatcre 
gradient is of the order of 50 K, with some uncertainty due to the fact 
that the single point of view of the pyrometer provides only an oblique 
view of the top and bottom. 

Sample outgassing can release large bursts of gas. In the work 
coil region, a pressure burst can raise the pressure to a level where an 
RF discharge can be supported. When this happens, safety interlocks 
must necessarily shut down power to protect the apparatus. Approximately 
two thirds of all fresh samples are lost prior to or near the melting 
point from this cause. In addition, a major portion of the lower tempera- 
ture deposition on the internal optics appears to occur from fresh 
samples. It is clearly necessary to preprocess all sami.>Ies before use 
to alleviate these problems. 

The calorimeter proper is a chrome-over-gold plated copper cylinder 
20.3 cm high and 13 cm in diameter. The center tapered receiving well 
is a removable copper sleeve 5-1/4 cm top inside dimension and 11 cm 
overall in length. A similar sleeve, wound with a manganin wire heating 
element, was used for electrical calibration. The calorimeter receiving 
well is lined with 0.010 inch tungsten foil and a variety of shards of 
tungsten stand in the well to prevent splash back. The splash protection 
has not been tested due to lack of liquid drops, but slightly surface 
melted samples have impacted the well with no damage and good heat 
transfer. This technique of splash protection has worked well for a 
variety of other systems [2;4j. 
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The total 25 beat capacity of the calorimeter as calibrated was 
7893 ± 7 J/^C. The basic uncertainty in total heat using this system, 
exclusive of sa^)le loss errors, radiant or electron beam power influx, 
but inclusive of reasonable jacket temperature excursions, thermometer 
errors, receiving well and foil liner weight errors, etc., is expected 
to be much less than 0.5 percent. 

A primary problem in this research effort has been the high tempera* 
ture measurement capability. Under 35 kV electron bomba rdsient, tungsten 
emits copious x*radiation. For this reason, visual access is highly 
restricted and a special imaging pyrometry system is needed. A General 
Electric silicon charge injection diode array camera is used to provide 
a 30 frame/sec oscilloscope ipiage of the sample and work coil area at an 
image magnification of about 5x. The levitated sample typically covers 
about half of the 10,000 diodes of the 100 x 100 camera area. A G£ 
constructed address/sample*and*hold unit allows selection of any one of 
the diode elements for analogue intensity output. The diodes have a 
dynamic range of about 20 (max. signal/avg. noise) and a linearity a. 
long term stab.'ity of better than one part in 10,000. The spectral 
response, however, (characteristically shown in figure 2) is quite 
non-linear. A major portion of this task effort has been devoted to 
implementing techniques necessary for conversion of this device into a 
pyrometr ically reliable tool. The initial effort by GE was to simply 
extrapolate the optical response of a selected diode by an exponential 
function from ribbon filament temperatures to the region of temperature 
interest. Given the fixed reference point of the melting point of 
tungsten, this technique should be good to 2 or 3 percent near the 
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reference point. As will becone apparent in the following discussion, 
this extrapolation is potentially poor, especially if the fixed reference 
point is not well identified. The probable error without considerable 
precaution could be of the order of 5 percent, i.e., perhaps as much as 
200 K. A reasonable goal in enthalpy determinations is one to two 
percent overall. This requires significantly better than one percent 
pyroflietric accuracy to achieve this error level in the face of emissivity 
uncertainties, sample temperature inhomogenieties , and the general 
problem of optics coating occuring due to sample outgassing and tungsten 
deposition. 

The optical path arrangement for pyrometry currently in use is shown 
ia Figure lb. Tfte consequences of restricted access, the camera charac- 
teristics, vapor deposition problems, and the lack of black body condi- 
tions to the implementation of accurate pyrometry require consideration 
of the theory and practice of optical pyrometry. 

Plank's law provides the basis for pyrometric measurements, relating 
spectral radiant flux to temperature by 

c^k~^/n 

~ exp (c^/ AT) 

where c^ and c^ are the first and second radiation constants (c^ = 
3.741832 X 10 W*ro; = 1.438786 x 10 ^ m*K), A is the spectral 

wavelength, and T is the absolute temperature. ^s then the radiant 
flux per unit area, ^r unit solid angle, per unit wavelength. The 
basis for the current International Practical Temperature Scale 
(IPTS-68) (16] at temperatures above the gold point, 1336.15 K is the 



ratio of the flux at teaperature^ T» to that from a body at some refer- 
ence te^>erature, T(ref). Figure 3, using the right scale, shows the 
temperature behavior of this ratio for a fixed reference temperature of 
1600 K and an ideal filter with unit transmission in a pass band 1 nm 
wide centered at the wavelengths given. This technique is usually 
referred to as brightness pyrosietry. 

Traditional pyrometry operates at a center wavelength of 650 nm 
with a typical band pass (including the observer's sensitivity) of about 
350 to 400 nm. Photoelectric pyrometers typically use band pass filters 
with a transmission window 35 nm wide. It is clear from Figure 3, 
that in the realm of 2500 to 4000 K temperature measurements, signifi- 
cently more sensitivity in temperature measurement is possible if shorter 
wavelengths are used. In this temperature region, there is still, 
however, more than a factor of two difference in absolute flux, favoring 
the longer wavelengths. At 3700 K, a 25 K change in temperature produces 
an absolute flux change of 6 percent at 450 nm, 5 percent at 550 nm and 
4 percent at 650 nm. 

At lower temperatures, another form of pyrometry has been used. 
Calleci two-color or ratio pyrometry, it uses the ratio of the fluxes at 
two wavelengths as the measured function of temperature. For any given 
wavelength pair, the sensitivity fails rapidly with increasing tempera- 
ture (see Figure 3, left scale and ratio curves). The high sensitivity 
' curves use wavelengths where the flux at the shorter wavelength is very 

j low at low temperature. However, a multicolor measurement using 

appropriate band passes can provide redundant temperature measurement. 
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Further, to the extent that eaiesivity on optics absorption is just 
independent of wavelength the measured temperature is independent of 
emissivity and changes in absorbance. With multicolor redundancy, it is 
possible to measure emittance directly. 

The general spectral response function of a silicon diode, shown in 
Figure 2, indicates a rapid fall off in sensitivity at shorter wavelengths, 
and the catastrophic termination of response at larger wavelengths, 
where silicon becoaies transparent. A choice of wavelengths is available 
in this region and it is not at all clear that the traditional value of 
650 nm is necessarily the best choice. 

The primary problems plaguing the spectral flux or brightness 
pyrometry method are the uncertainty in emissivity and absorption effects. 
The technique obviously requires complete knowledge about the actual 
transmission of all optical train elements, but is relatively less 
sensitive to the spectral behavior of the absorbances in the path as 
long as they do not change and the pyrometer* s spectral response band 
pass is sufficiently restricted. 

In order to implement the brightness pyrometry technique, it was 
necessary to provide a reference source. Traditional pyrometers super- 
impose an ioiage of an internal standard lamp on the attenuated image to 
be measured. The standard lamp flux is then adjusted to match the 
attenuated source image and the ratio is accomplished as a null measure- 
ment. The range of comparison between internal standard lamp and external 
source generally uses a broad band absorber as the attenuator. The 
external comparison of a black-body sorree and a hotter body is accom- 
plished by attenuating the light from the hotter body by a device called 
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a sector wheel; essentially this is a chopping wheel behaving as a 
Mechanical aperture utilizing the eye as a visual snoothing device. 

This Method can provide for a direct reading instrument of excellent 
accuracy and stability, but has no capability for dealing with changing 
absorbance in the optical path, such as occurs when vacuum deposition 
occurs. Vacuum deposition on the vacuum pyrometer window and reflection 
mirror changes the transmission of the optical train typically by a 
factor of two to five in the course of a single experiment. This magni- 
tude of change is equivalent to a temperature change of 300 to 500 K. 

The installation of a reference lamp internally within the vacuum 

system was thus required. The mirror. Ml, shown in Figure lb, reflects 

an image of the sample through the vacuum window which is located so 

that no direct line of sight exists between the uample and that window. 

This essentially eliminates coating of the window and provides x-radiation 

protection. The small mirror M2 superimposes an image of the internal 

reference filament at the image plane of the sample, offset to one siae 

about 1/2 cm visually from the sample limb. M2 shades this filament 

from deposition and its orientation shields its own surface. Mirrors Ml 

2 

and M2 are 1/4 wave flat first surface mirrors, 7 cm as used, obtained 
from Edmund Scientific Corp. The installation of a General Electric 
18A/6V/T10 strip lamp provides the means to establish an initial tempera- 
ture scale, by comparison with a Pyro micro optical pyrometer, and the 
means to make a running measurement of deposition on mirror Ml. 

The p cedurc which has been developed involves replacing the 
levitation coil assembly with a thicp glass plate and pumping the system 
down when a new lamp must be installed. The lamp is then operated over 
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the region from 1600 to 2450 K, brightness temperature 1750 to 
2770 K true tesiperature) , while measuring the current drawn by the lamp 
as a function of brightness temperature monitored by the pyro optical 
pyrometer observing through the glass plate. The current is measured 
with a 4^ digit voltmeter monitoring the voltage drop across a high 
current shunt. Current is supplied to the lamp from a well regulated dc 
supply. At the same time, the GE camera also monitors the strip lamp 
through its optical train. The emissivity of tungsten selected [17] was 
0.425 for this filament. The procedure must be repeated frequently, as 
the lamp reference is operated far above its long term stability point [18] 
This technique provides a low temperature reference to the IPTS-68 once 
the absorbance of the thick optical plate is corrected. This is accom- 
plished in a separate experimt by measuring the temperature change, 
as a function of lamp brightness temperature, with and without the glass 
plate in the optical path. The sighting path through the plate was 
identified by physically markirig the plate and carefully measuring and 
preserving the angles of the calibration experiment. The preferred 
method of performing this correction is to actually measure the 
absorbance [4] but in most cases, an actual temperature effect 
measurement is adequate. The Wien approximation to the Planck Law 


c^\^/n 

^ expCc^/XT) 


is a good approximation in this temperature wavelength regime. Since 


transmittance is just 


I -ax 

•j- = e = transmittance 

0 
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and in the Vein approximation, behaves just as an emittance, identifying 
the flux change as a transmittance effect gives the re.lation between 

brightness temperature Tg observed without the absorber to the observed 

•k 

temperature T ^ with the absorber as 

apparent 


k - 1 

k k 
T T 

B apparent 


= ^ ln(e ^ 


absorber 


The glass plate used showed a variation in K of a factor of two over the 
temperature range 2300 to 2500 K. This is attributed to a color depen- 
dence in the glass transmission, and the values were applied 

at the Tg values measured. The effect of this absorber was to lower 
observed brightness temperature by 95 to 52 K. The error due to uncer- 
tainties in the exact wavelength dependence were estimated from the 
direct lamp temperature versus current measurements to be of the order 
of 15 K. This plate must be replaced in future measurements by a sheet 
of pyrex or quartz whose visible spectral characteristics are known and 
well behaved. 

The imafng array camera uses a relay lens with an effective aperture 
of 1 in. to provide the necessary magnification to nearly fill the frame 
with the sample image. To convert this device to a pyrometric device, 
it was necessary to control in a known way, both the relative sensitivity 
and the true spectral wavelength response of the system. As a result of 
early trials, a combination filter holder, aperture control device was 
designed and put into service. This device fits on the entrance aperture 
of the relay lens and has two manually turned detented wheels. One 
wheel holds one inch aperture narrow band 80 mn FWHM) interference 
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filters blocked to 10 or better outside the baud pass. The two filters 
actually calibrated in this work had band passes centered at 651 nm and 
451 nm. An auxiliary broad band pass neutral density (ND) filter mounted 
on front of the filter/aperture device was used as a range control. The 
other wheel has a series of knife edged circular orifices: 1, 5/8, 1/2, 
5/16, 3/16, and 1/8 inch in diameter, which were to be prepared to 
± 0.0005 tolerance. 

Calibration consisted of comparing the output of a specific diode 
of the imaging array as a function of temperature, filter aperture and 
ND filter with clean freshly mounted mirrors using the internal cali- 
brated lamp. A variety of measurements for cross checking redundance 
were taken to allow the calculation of any observation from others. The 
one inch aperture was found to be larger than the system aperture, as 
expected. The 5/16 aperture consistently showed itself to be 7 percent 
too small in area or 0.010 inch smaller than the manufactured 5/16 
(0.312 inch) aperture. On subsequent check, the actual diameter was 
found to be 0.302 ± 0.001 inch, exactly as predicted. All other apertures, 
and the transmission factors of the two filters agreed with actual 
temperature measurements with a reproducibility of better than 0.6 percent. 
The quality of the fit of data over the lower temperature calibration 
range fit the Planck Law expression with a maximum error of 15 K and an 
average error of 5 K. When translated by aperture changes to the tempera- 
ture range 3500 to 4100 K, cumulative errors in aperture and fixed 
optical path absorbances are expected to contribute less than an additional 
10 K error. With the built-in reference, the remaining problem is 
changes in the spectral reflectivity of mirror Ml. In fact, the very 


77 


effort to measure the problem, taking pyrometric measurements at two or 
more wavelengths, allows an independent check on temperature when the 
wavelength dependence on reflectivity of Ml is measured. Where possible, 
experiments are arranged to attempt this, but the inability of the GE 
system to provide a significant time of constant temperature liquid 
levitation (or even hold on to the completely liquid sample) has restricted 
these measurements to stably levitated solids. In that case, the agreement 
between the three temperatures, brightness at 650, brightness at 450, 
and color ratio between these two measurements has a precision of better 
than ± 20 K, indicating accuracy attained of 0.6 percent. 

DISCUSSION AND CONCLUSIONS : 

The main goal of this phase of this task has been to help RICE install 
a calorimetric subsystem on GE^s existing and proven [19] liquid tungsten 
levitation system. The necessity of converting an approximate temperature 
tool into a precision pyrometer was anticipated, and although a variety 
of problems were encountered (such as the need to operate around extreme 
radio frequency interference (RFI) from adjacent unrelated experiments, 
internal problems with the diode select logic, etc.), the procedures 
which have been established are workable and meet the initially antici- 
pated needs of the project. The calorimetric system has performed 
perfectly. The addition of a data logging system at the bepiuning of 
this phase was made eliminate the need for an observer for the calorim- 
eter heat trend. As expected this addition allowed thr expansion of 
effort from two to as many as six drop attempts per day. The 8 kJ heat 
capacity of the calorimeter block and the vacuum insulation still restricts 
the effort to one sample in the calorimeter in four hours. 
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It had t.^en ar^iclpated that if this phase progressed as in fact it 
did, that we would have obtained several successful liquid drops and be 

in a position to refine the current value uncertainty in the heat of 

I 

\ fusion, H (tungsten). 

I m 

I The Tamman Rule [5] (AH /T = S -v 2.3 cal/mol-K) vaV. e of 8.5 

t m m 01 

f 

I kcal/mol has been the operational value in all the standard references 

[9;20]. Dikhter and Lebedev [21] have obtained a value of by an 

exploding wire technique. Their direct value, 80 ± 4 cal/g, corresponds 
to a value of 14.7 ± 0.7 kcal/mol. This value is remarkable as is 
implies a value of AS =4. The largest value of AS known for the 
transition metals is molybdenum, with AS^. = 3.08 cal/mol*K [3a]. Even 
though the general trend in AS^ seems to be increasing down each group [4] 
and tungsten is below molybdenum, this value still appears to be too high. 
Kubachewski, et al. [22] give a value of 12.1 kcal/mol for All^(silicon) 
which was expected [23] to have the * ighest heat of fusion of any element. 
Vouch, et al. [19] report a value for AH^ of 11 ± 1 kcal/mol obtained by 
applying a cooling curve model to experimental data. The model ignores 
supercooling, but provides a value of AS^ of 2.98 ± 0.3 cal/mole*K, in 
line with expectations from the molybdenum data. Tne value is still 
very uncertain, as Vouch, et al. [19] repeatedly noted, especially in the 
weakness of this derivation for aon^spherical samples. (Terrestrial 
levitation). In the course cf ^his work, to establish the HE ability to 
levitate and melt, a sample was allowed to fall onto a copper plate in a 
one meter drop tube. The sample splashed as though completely molten, 
and all (± 1 percent) fragments were recovered. Since the temperature 

was estimated to be just above melting, even if AH^ is 8.5 kcal/mole, 
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significant supercooling must have occurred and the Wouch, et al. (19] 
model for the cooling curve type of determination [24] should not apply. 

Unfortunately* even though the calorimcitry and pyrcmetry are now 
functional, the levitation process itself has been plagued with problems. 
Although outside the area of this task, improvement in this area is 
essential to the task and efforts to deal with the problems have ^een 
made. The work coil was changed at our suggestion. The problem of the 
so-called "hot short", where samples would deform, "melv", or explode at 
temperatures hundreds of degrees below the melting point of tungsten was 
again encountered. Previously, this problem had bec*i diagnosed as being 
caused by hydrocarbon contamination durir^ machining of the specimens. 

This recurrence was identified as vapor deposition of volatile metals 
from the soft solders used in the coil area. Figure 4 shows the effect 
of this problem on a group of samples. In addition to the loss of 
samples, the shape change can be mistakenly interpreted as melting, 
making the use of visual observation of "melting" as an internal p>ro- 
metric calibration point unreliable. 

The most serious problem has been the inability of the levitation 
system to melt and hold tungseen. More than 30 samples have been luvi* 
tated and successfully e-beam heated in full dress calorimetric runf 
In no run has it been possible to contain the liquid for more ther 1 to 
2 seconds before inadvertant pouixng from the coil occured. With the 
current level of hand operations necessary, the very short, unpredic- 
table containment times give an intolerable failure rate. The consequences 
of the levitation containment problem are quite importan\ the eventuil 



success of this progras. The current inplication under consideration is 
that the sirface tension/density ratio of liquid tungsten is perhaps as 
low as that for gold or lower. This problea was encountered at RICE in 
an effort to neasure enthalpy increments for tantalum [4], The measure- 
ment effort for platinum, also required special efforts to contain. It 
is quxte possible that this will he a general property of the entire 
third long period of the transition elementr. Improving containment 
could require a significant additional effort on the part of GE to 
improve their levitation facility 

Because of the containment pre , it has been lecessary to suspend 
experimentation in this effort. At a minimuir, ^t vill be necessary to 
automate the calorimeL:>.ic system to ijiprove the system response time to 
the point where the ability of the calorimeter gate structure to respond 
no longer requires long or predictable containiLv:nt times. The general 
characteristics of such a system have already been examined and will be 
the subject of future reports when the various aspects have been discussed 
fully with RICE and GE. RICE is expecting a new post-doctoral student 
to enter this 7 ^oject and the final form of the plans will, of course, 
depend on this. 

The advantages of an automated calorimeter system extend beyond the 
inErediate application. It is expected that calorimetric studie.r will be 
among those sele -ed for roicrogravity experiments. The current effort 
requir«_s thiee or *^cur researchers tc operate the a^oaratus and take 
data. The reduction of dat? ueat content and temperature requires 
some rime as well. The effort of automating the calorimetric systep and 
Integrating the pyre aetry system into such automation *tviously is a 
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step which will be done in the no^aal progress of this experimentation. 
This seems to be the obvious time to implement such automation. It is 
recommended in the interest of future efforts, that the pyrometric 
system be automated as well. The major item necessary to providing an 
automated pyrometric system is some method of measuring the change in 
absorbance of the primary mirror. Ml (Figure lb). 

Figure 5 shows a schematic of the hardware portion of a proposed 
automatic imaging pyrometric system. The motor driven aperture and 
filter assembly on the pyrometer is a simple extension for computer 
control. The novel feature is the modulated photodiode assembly. With 
a chopping wheel modulating the internal filament, a simple photodiode 
ccn be used to extract the change in apparent brightness of the reference 
filament by lock-in aaq>lifier techniques. With this information in 
real-tiiie as wel^ as the output from the imaging array, a microcomputer 
can, with careful calibration as has already been developed and described 
in this report, provide an almost real-time display of true temperature, 
and report the onset of problems in temperature measurement from the 
multiple redundancy inherent in mulcicolor pyrometry. 

The same computer can easily have the capability to handle the calo- 
rimeter autofliation and data logging requirement. The same data system, 
with auxiliary floppy disk storage, can provide complete data reduction 
in a few minutes after the experiment ends. In fart, a running luonitor 
of the final period cf the #*xperiment may be able to significantly 
shorten the 4 hour p^r run time period currently necessary for best 
work. 
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FIGURE 1q 


Fig. la. Front cross sectional V3iw of the C£-Rice roupied apparatus 
showing the relative orientation of najor components. 
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FIGURE 1b 
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A top view showing details of the optical arrangement for pyrometry. 
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Fig. 2. Characteristic spectra! response behavior for silicon photodiode 

type devices. The GE imaging array elements have a similar response 
function. 
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SPECTRAL WAVELENGTH RATIO 


PYROMETRY TEMPERATURE FUNCTIONS 



TEMPERATURE (K) 


Fig. 3. Plank law functional behavior for two type of pyrometry. The 
dashed (and solid) curves are spectral radiance ratios ('right 
scale) for the labeled wavelengths for = 1600 K, typical of 

brightness pyroiretry ♦ The dot-dash curves show the flux ratios 
(left scale) for wavelength pairs, typical for multicolor pyrometry. 
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Schematic of opt i ca 1-mechanica 1 arrangement for a computer operated 
real time direct reading pyronietric system. See text for discussion. 
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High Temperature Thermophyslcal Properties of Refractory 
Materials from Free Cooling Experiments 

Jack H. Colwell, Ared Cezalrllyan, and Lawrence A. Schmid 

Thermophyslcs Division 
Center for Chemical Physics 

Summary 

Free cooling experiments are to be used to obtain high temperature 
(1200-4000 K) , heat capacities thermal radiative properties, and thermal 
dlffuslvltles of materials while samples are levitated in a space environment. 
The data are to be derived from the simultaneous measurement of the total 
ladlant hewt loss and the temperature change of spherical samples as they 
freely cool by radiation alone. Pyroelectric detectors are being investigated 
for use In making the total radiance measurement. The dynamic response of 
these detectors have been tested using a millisecond pulse heating apparatus, 
which permits metal straps to be heated at rates up to 4000 K/s and cooled 
at rates up to 1500 K/s. In these experiments the pyroelectric detector 
measures the total normal radiance from an exposed portion of the strip 
while the radiance temperature of the strip 1& determined by a high speed 
pyrometer. A small (^1%) discrepancy between the heating and cooling 
portions of the experiments Is ascribed to the radiance of the strip, not 
to the behavior of the detector. Normal total emittances derived from 
these hl^h speed measurements are in excellent agreement with t^e steady 
state results of others. The large temperature gradients that can cevelop 
In rapidly-cooling spherical samples complicates the Interpretation of the 
free cooling experiments. Mathematical expressions have been derived for 
values of the physical parameters of materials from the observed time 
depended e of a sample's surface temperature and total radiance. 


91 



I 


I 

I 

Introduction 

The objective of this research is to measure the high temperature (1200^ 
4000 K) thermophysical properties, viz., heat capacities, th^^rmal emittances, 
and possibly thermal dlffusivltles, of refractory materials. The experiments 
will exploit the mlcrogravlty environment of space to maintain the samples, 
when hot, without physical contact. RF induction will be used to heat samples, 
and subsequently, will be used at low power levels to maintain the position 
of the samples. Samples will be spherical, fabricated in that form as solids 
and, when llqui:^, maintal**ed in that form by their surface tension. Experiments 
will consist of heating the samples to a high temperature in a cold-walled 
vacuum chamber and then allowing the samples to freely cool by radiation. As 
the sample cools, its surface temperature will be determined by multicolor 
radiation pyrometry and the total radiative heat loss will be determined 
simultaneously by a wide-band pyroelectric detector. The spherical geometry 
of the sample makes it possible for the total radiance measurement to be made 
with a single detector. A unique feature of these experiments is that it 
will permit liquid samples to be followed into the supercooled region all the 
way to their point of recalesccnce, thus yielding thennophysical data on 
supercooled liquids and possibly kinetic data on the solidification process. 
Containerless liquid refractory metals undergoing free cooling may be ex, :,ted 
to supercool by as much as 1000 K. 

The cooling rates of the samples will be a functi^ of the temperature 
and the size of samplest ^ » 2] ^ yov smaller samples, less than 5 mm diameter, 
at 4000 K, the cooling rate will exceed 1000 K/s. Larger samples, 10 mm 
diameter at 1500 K will be cooling at rates of the order of 10 K/s. 
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Because of the dynamic nature of the cooling experiments, temperature gradients 
will exist within the samples [^1. Frr metals with a high thermal diffusive ty, 

**he total ten 4 >erature differential from surface to center will be a maximum 
of tens-of-degrees and will have little effect on the measurements. For 
materials of lower dlffuslvlty, the temperature differentials can reach 
hundreds-of -degrees. The temperature differential is, however, an Inverse 
function of the sample radius, so by measuring different sized samples the 
physical parameters can still be determined. For materials within a particular 
range of dlffuslvltles, the dlffuslvlty itself can be obtaxnea. The mathematical 
procedures for determining the values of the material parameters from the 
experimental observables is described in detail in Part B of this report. 

One of the innovative aspects of this work is the use of pyroelectric 
detectors in dynamic measurements of the total radiance, in Part A of this 
report we describe the experiments with a pyroelectric detector using the 
rapid electrical heating apparatus [^1 of the Dynamic Measurements Laboratory 
at NBS. In these experiments, metal strips are heated to a glvei> temperature 
in less tu ■ one second, and then allowed to cool. The temperature is 
measured pyrometrlcally and a portion of one surface is viewed by the pyroelectric 
detector. The pyroelectric detector gives a determination of the total 
normal emlttance of the surface as a function of temperature. The main 
advantage of this apparatus is that heating and cooling rates can be generated 
which are several times faster than expected in the free-cooling experiments, 
thus thoroughly testing the dynamic response of the detector. 

The vacuum assembly is nearly complete and the inductive power supply 
is now on hand for the ground-based simulation of the free-coolinp experiments. 
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This effort is being detained while we wait for the renovation of a new 
laboratory which will house the experiments. In addition to the experiment 
described in this report, ve are proceeding with the development of a 
millisecond- resolution two-color pyrometer which will be tested in the 
simulation experiments. 



Part A: 


Experimental Procedures 

The millisecond pulse heating apparatus^ used in these experiments 
in^resses a nearly constant voltage across a conductive specimen for a 
predetermined length of time. The constant voltages are produced by a bank 
of batteries. Voltage levels are determined by the number of cells in the 
circuit and an adjustable series resister. By the time the sample reaches 
the temperature where measurements can begin, T > 1200 K, all switching 
transients have died away and the current is constant except for the slow 
decrease caused by the increase in the resistance of the specimen as its 
temperature rises. The duration of the pulse is determined by the circuit 
characteristics and a temperature limit, and generally is less than one 
second. 

The sample was viewed from one side with a high speed pyrometer which 
fo::used on a target area of 0.2 mm diameter at the center of the specimen. 

The pyrometer uses a photomultiplier tube as a detector, which is alternately 
exposed to light from the specimen and from a calibrated strip lamp. The 
sequence of measurements made by the pyrometer is govemeJ by a 12 aperture 
chopper rotating at 200 Hz. Six apertures admit light from the specimen and 
are completely open as are two of the apertures admitting light from the 
lamp. The other four apertures contain two pairs of light attenuators. With 
each rotation of the chopper, six observations of the specimen are made, 
interspersed vith six reference observations of the lamp. The pairs of 
reference attenuators used in the cho^ r are not perfectly matched so their 
signals show up as two closely spaced doublets on the output (see Fig. 1). 

The specimen temperature is determined only when the signal lies between 
the upper and lower reference lines; this constitutes a change of 20^25 
percent in the brightness temperature of the specimen depending on 
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the particular temperature range. Temperatures in different ranges are 
obtained by placing calibrated neutral-density filters in the lieht path 
from either the lamp or the specimen. To cover a vide range of temperature, 
therefore, successive experiments must be run with different filter combinations 
to cover each successive temperature range. 

The electrode assembly that holds the specimen is enclosed in a cylindrical 
vacuum chamber that is sealed by an 0-ring to the bottom plate. Thr pyrometer 
views the specimen through a glass window in the side of the vacuum chamber, 
with the objective lense 15 cm from the specimen. The pyroelectric detector 
is on the opposite side of the specimen, mounted in a side tube of the 
vacuum chamber so that no windows need intervene between the specimen and 
the detector. The light chopper, which is required for pyroelectric detection, 
is a small (30mm diameter) rotating sector device and is operated at a 
chopping frequency of 200 Hz. The chopper motor has teflon bearings 

and has now been operated for tens of hours in a vacuum with no apparent 
overheating or other ill effects. 

The portion of the specimen viewed by the pyroelectric detector is 
determined by a pair of horizontal knife edges that obscure the top and 
bottom portions of the strip, so that the detector views only the central 
portion. The knife edges are mounted on the electrode support assembly, 5 
cm from the specimen. The detector is 23 cm from the specimen so the 
effective sample length is 25% larger than the slit width. The active area 
of the detector is 1 mm in diameter so the penumbra region increases the 
effective length of the specimen by 0.2 mm. In these experiments, the gaps 
between the knife edges were 12 to 15 mm, so any error in the penumbra 
correction should have a negligible effect on the measurement. Our first 
experiments were with tungsten strips, 6.35 mm wide, 0.25 mm thick, with 
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45 on of length exposed between the electrodes. The opening between the 
knife edges was 14.6 nn. Tlie outputs of the pyrometer and pyroelectric 
detector were recorded with two dual-trace oscilloscope-camera combinations. 
Numerical values were obtained from the photographs using a traveling microscope. 
A subsequent set of measurements was made with niobium strips of the same 
size as the tungsten strips and with a knife edge opening of 12.7 mm. For 
this second set of measurements, data were recorded us'^rg a pair of digital 
oscilloscopes from which the data could be transferred directly to the PDP-11 
computer for processing. The light chopper for the pyroelectric detector was 
synchronized with the pyrometer chopper so that the two sets of date could be 
accurately correlated. 

The linearity of the response of the pyroelectric detector was established 
using the Inverse distance law. A DXU 1000 W tungsten halogen lamp was used 
as a light source. The detector was moved over distances from 50.0 to 
141.4 cm from the lamp. Measurements were made for Irradlances covering 
more than three orders of magnitude from 32 mW/cm^ to 20 pW/cm^. After values 
were measured ^rom 50 cm out to 141.4 cm, the lamp power was reduced until 
the detector signal was the same at 50 cm as it was previously at 141.4 cm. 

The process was thei\ repeated. These tests showed that the linearity 
of the detector was better than '% over this wide range. The principal 
uncertainty in the caLlbratlon procedure comes from the geometric 
distribution of the filament in the light source. The filament of the lamp 
is a colled-coll construction, approximately 20 mm long and 5 cm in 
diameter. There is some uncertainty in the point taken as the origin 
for the distance measurement. Also, the angle of the lamp to the light 
path to the detector has an effect, probably due to changes in the amount 
that the front colls of the lamp obscure tl\e back colls at different detector 
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distances. Care has to be exercised ::hat the lamp had reached thermal 
equilibrium after each power change. 

Experimental Results 

Oscillographic traces of one of the tungsten strip heating experiments 
are given in Fig. 1. The two oscillograms are of the same experiment with 
the data collected at different sweep rates on separate oscilloscopes. The 
upper trace in each oscillogram is the pyrometric temperature signal, 
increasing to the point where the current is cut off and then decreasing as 
the specimen cools. The horizontal bands are the pyrometer reference 
levels, their corresponding radiance temperatures are indicated on the 
lower oscillogram. The lower trace on each oscillogram is the AC signal 
from the pyroelectric detector, the peak-to-peak distance being a measure 
of the total radiant energy. The shape of the pyroelectric signal was 
puzzling when first observed. With slowly varying light intensities, the 
two halves of the AC signal are found to be symmetrical, so the asymmetric 
signal observed in the rapid heating experiments was unexpected. The upper 
half of the pyroelectric signal in the oscillograms corresponds to light 
falling on the detector, and the lower half to the detector obscured by the 
light chopper. A clue to the explanation of the asymnetric signal is found 
in Fig. 2 which is the response of the detectoi when suddenly exposed to a 
constant radiation source. The initial response is a large positive signal 
during the exposed half of the cycle, and a very small negative signal 
during the obscured half of the cycle. With time the positive peaks decay 
and the negative peaks increase until they are symmetrical about the zero 
signal axis; the relaxation time, t, for this process is approximately 100 
ms for this detector, r is the thermal relaxation time of the detector and 
is determined by its heat capacity, c, and the thermal conductance, h. 
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between the detectc'r element and Its mounting, so that t “ c/h. The Initial 
response In Fig. 2 la when the detector la first heated above asdilent, the 
temperature Increase la large In the exposed cycle and small In the obscured 
cycle because little heat Is being lost to the heat sink. As time goes on, 
the average temperature of the detector Increases and more heat Is being lost 
to the heat sink, so that the temperature Increase during the exposed cycle 
becomes less and the decrease during the obscured cycle becomes greater, 
eventually becoming equal when the steady state Is reached. For a steady 
light source, as In Fig. 2, this relaxation phenomena Introduces no detectable 
error (~0.3Z resolutlor In this experiment) between the peak to peak voltages 
at the beginning of the relaxation period and those after the steady state 
was reached. 

A second relaxation phenomena that affects the signal from a pyre electric 
detector Is the rise time of the device. For the detector used In these 
experiments the time constant was approximately 0.5 ms. As a consequence of 
the rise time and the drop In the signal due to the thermal relaxation, the 
maximum signal with constant Irradiation occurs at approximately 2.0 ms Into 
each half cycle. 

The series of experiments on tungsten strips covered thermodynamic 
temperature ranges from 1500 to 3000 K. Heating rates varied from 4000 K/s 
at the lowest temperatures to 2000 K/s at the highest temperatures and the 
corresponding cooling rates ranged from 500 uo 1500 K/s. A direct comparison 
of the normal total radiance with the radiance temperature of the specimen 
showed that the measured total radiance when the specimen was cooling was 
always 1 to 2X lower than when It was being heated, except for the region 
near the maximum temperature. 
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Themodynanlc tenperatures were calculated from the measured radiance 
ten 4 >erature 8 using the spectral emittance data of de Vos^^^ and Abbott^ 

The Stefan*Boltzmann law and the total normal radiance data were then used to 
calculate the total normal emittance of the tuiigsten strips as a function of 
temperature. The calculated rvalues agree within 1% with the publish values 
of Abbott^ which were obtained In steady-state experiments. This close 
agreement on an absolute scale Is fortuitous for the Inte.itity calibration of 
our detector was not tho'ight to be that precise. 

Fig. 4, and 5 are the results of experiments with niobium strips over 
three different temperature ranges; these are direct computer plots with no 
Intermediate treatment of the data. The plot of normal total radiance 
(pyroelectric voltage) versus radiance temperature again shows the slight 
difference between the heating and cooling curves. The difference increases 
the longer the sample cools and is largest in the higher temperature ranges. 
Normal total emittance values from these experiments are within 5-10% of the 
steady-state values determined by Abbott the discrepancies resulting from 
the lack of precise values for the normal spectral emittance required to 
derive chermodynamic temperatures. 

Fig. 6 presents the normal total radiance and radiance temperature 
measurements during the melting of a niobium strip. The radiance temperature 
shows the flat plateau, characteristic of melting cur^'es in these dynamic 
experlmentsf^l ; the radiance temperature drops by only 2-3 K in a period of 
200 ms. The normal total radiance, however, has a short plateau and then 
drops off steadily (down 3-4% after 100 ms) until it finally drops precipitously 
as tne specimen melts through. In an earlier melting experiment the timer 
switch stopped the current prior to the specimen melting through, the melting 
plateau was only about 100 ms wide. The measured width .rf this partially 
melted sample after the experiment was found to be 1-2% narrower 
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than It was originally btcauae of the melting at the edges. This 
Indicates that the fall off In the normal total radiance In Fig. 6 may 
be due to a narrowing of the specimen. The radiance temperature Is 
unaffected by the narrowing of the apeclmen because It Is determined by 
observing a small spot of the center of the strip. 

The melting point of niobium Is known, 2750 so the normal 

total emlttance of liquid niobium at the melting point can be evaluated; 

It Is 0.256. 

Conclusions and Discussions 

The small discrepancy between tne heating and cooling portions of 
the correlation between the pyroelectric voltage and the radiance temperature 
is thought to due to changes In the total radiance determinations and 
not due to any error In the temperature measurement. We feel that the 
change is most likely due to differences In the total radiance of the 
strip during heating and cooling. While cooling, the strip will be 
losing energy from the edges of the strip (4% of surface area) as well 
as the faces so that the average radiance from a face will be less than 
the radiance from the center of the face. During the heating period 
this tendency is mitigated in two ways. The rate of heating is much 
faster than the rate of cooling, so that the temperature gradients will 
not be as fully developed. More Importantly, however, when heating, the 
current will tend to follow the path of lowest resistance, l.e. the 
coolest path, and will suppress the development of the temperature 
gradients. There is the possibility that the discrepancy is inherent in 
the pyroelectric detector. Because of the finite rise time (and the 
thermal relaxation) of the device, a different voltage signal could be 

obtained for rapidly increasing and decreasing signals. This possibility 
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has not been Investigated, but If it were a najor contributing factor to 
observed discrepancy. It would be expected to be largest near the 
of the teaperature curve were the signals are changing most rapidly; 
this £& not the case. In any case, the heating and cooling rates in 
these strip heating experiments are several times larger than those that 
will be encountered In the free-'cuollng experiments, so errors from this 
source will be minimal. 

The principal aim of this stage of our research Is to establish that 
pyroelectric detectors can be used to make accurate total Irradl^nce 
measurements In high speed experiments. These experiments have accomplished 
that. 
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Fig. 1 Tungsten strip heating experlneat. Two osclllograiiis ol the 

sane experlisent at different sweep rates, SO ns/cn in the cop 
oaclllogrsn, 200 ms /cm In the botton. The top curve In each casi'i 
Is the pyrometer output, the bottom pair of curves is the ac 
pyroelectric voltage. The horizontal traces are the reference 
levels of the pyrometer and their corresponding radiance temperatures 
are Indicated on the lower oscillogram. The radiance temperatuT e at 
the peak is 1912 K which corresponds to a thervodjrnamlc teisperature 
of 2056 K. 
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Fig. 2 The pyroelectric response upon Instantaneous exposure to a 

constant light source. Sweep rate Is SO ms/cm, pyroele trie 
chopping frequency Is 200 Hz. Variation of the two voltage levels 
is due to thermal time constimc of the detector; the peak-to-peak 
voltage is not affected by the relaxation process. 
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Niobium-Strip heating experiment to 1900 K. Arrowy o;i curv s 
indicate scales to be used. The pyroelectric voltage (proportional 
to normal total radiance) and the pyrometrically determined radiar^’C 
temperature are shown as a function of time, and also plotted 
against each other. On the latter curve the encircled points are 
those with the temperature increasing. Near the maximum temperature, 
the heating rate is 4200 K/s and the cooling rate is 500 K/s. 
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Fig. 4 Niobium-strip heating experiment to iflOO K. Arrows on curves 

indicate scales to be used. The pyroelectric voltage (propc>rtional 
to normal total radiance' jind the pyrometrlcally determined radiance 
temperature are shown as a function of time, and also plotted 
against each other. On the latter curve the ep-'ircled points are 
those with the temperature increasing. Near the maximum temperature, 
the heating rate is 3600 K/s and the cooling rate is 700 K/s, 
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Fig, 5 Niobium-Strip heating experiment to 2375 K. Arrows on curves 

Indicatescales to be used. The pyroelectric voltage (proportional 
normal tota radiance) and the pyrometrically dotcrmined radiance 
temperature are shown as a function of t5jne, and alsc plotted 
against each other. On the latter curve the encircled points are 
those with the temperature increasing. Near the maximum temperature, 
the heating rate is 200 K/s and the cooling rate is 1200 K/s. 
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Fig. 6 Nlobitiai-strlp melting experiment. The pyroelectric voltage 

(proportional to normal total radiance) and the pyrometrlcally 
determined r&dlance temperature are shown as a function of time. 

The plateaus of the curves Indicate the period where the surface of 
the specimen Is melting. The precipitous drop In the pvroelectric 
voltage indicates the point that the specimen melts through and 
collapses. The melting point of niobium Is 2750 K. 
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Part B: 


Mathenatlcal Analysis of the Free-Cooling Sphere 
1. Overview of the Axxalysls 

The free-cooling calorAmetrlc method discussed In this report alms at the 
determination of the temperature dependence of the thermal parameters (heat 
capacitance, dlffuslvlty, and emlttance) of refractory materials at high 
temperatures by optical observations of the time dependence of the surface 

temperature T (t) and the radiance E(t) (radiated heat flux per iinlt area) of a 

s 

hot solid sphere of the material as It cools down. (The possibility of melting 
will not be discussed In what follows. This will be the subject of future 
research.) 

The eml .tance e(Tg) of the surface Is simply given by 

e - E/(oT “♦) (1) 

s 

where o is the Stefan-Boltzmann constant. Thus the real problem Is to express 

the dlffuslvlty a and the heat capacity per unit volume c in terms of S(t) and 

Tg(t). In the limit a -*• “ (or R 0 where R Is the radius of the spherical 

sample) , the temperature throughout the sample can be taken equal to T , and c 

s 

can be found from the following simple relation: 

-(4TrR3/3)cdT /dt * 4 ttR2e or c = -3E/RT (2a, b) 

s s 

where T = dT /dt. For finite a, however, the temperature Inside the sphere is 
s s 

higher than T^, and so eq. 2b can only be the leading term In a more compli- 
cated expression concalnlng higher-order correction terms involving a. Such an 
expression for c will be derived in what follows, first for the case of constant 
c and a (cf. eqs. 30 and 37), and then for the case of temperature-dependent 
parameters c(T) and a(T) (cf. eqs. 56 and 57). Similar expressions will be 

derived for a as a function of the time derivatives of T and E (cf. eqs. 35 
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and 55). 


An essential feature of the free-cooling method described in this report 

is that the temperature field within the sample is assumed to be unknown at the 

beginning of the experiment. Thus the heater, be it an induction coil, furnace, 

cr some other device, need not be precisely controlled in order to assure that 

a specified quantity of energy be injected into the sample, or that a specified 

initial Interior temperature field be achieved. The price that must be paid 

for this freedom is that the readings of T^(t) and E(t) for a time interval 

At-R^/a at the beginning of an observational run after the heater is turned off 

cannot be used to determine a(T) and c(T) , because T^ and E are too sensitive 

to the details of the unknown Initial interior temperature field. It is only 

when the initial transients have died out and the thermal behavior of the 

sample enters the "post- transient regime" that tne expressions for a and c in 

terms of E(t) and T (t) (and their time derivatives) that are derived in this 
s 

report can be applied. For the sample sizes contemplated, however, only a 
couple of seconds are required for the transients to die out in the case of a 
poor thermal conductor like uranium dioxide, and this is to be compared with a 
toL<il observation time of nearly 30 seconds • Tn the case of metal sample? th#> 
time required for the transients to die out is a much smaller fraction of the 
total observation time. These estimates are discussed at the end of Section 3. 

In order to relate the analysis of this report to the existing ?.iterature 
on the analysis of free-cooling, it should be noted that there are three 
different kinds of cooling problem which must be clearly distinguished: 

(1) the predictive or initial-value problem in which the Initial temperature 
field within the body is given, and it is necessary to calculate Tg(t) 
knowing the functional dependence of the radiance on une surface 


] 10 


temperature; 



(1) nhc InvTift problea, In which E(T^) ia known, T^(t) la given, and it la 
lacaaaary to calculate the initial interior teaveratore field; 

<3, the caloriaetric problea, in which E(t) and T^(t) are given, and it la 

neceaaary to calculate c(T) and o(T). A knowledge of the initial interior 
teaperature field la neither given nor dealred. (Thla category of problea 
ia aoaetlaea alao referred to in the literature aa the "inverae problea".) 
This report will concern itaelf only with the flrat and third categorlea. Only 
the weak fora" of the predictive problea will be dlacuaaed, by which ia aeant 
the prediction of T^(t) (and the interior teaperature field) in the poat- 
tranaient reglae. Inaaauch as T^(t) in thia raglM ia deteradned by the theraal 
pnriiaatera alone, it ia poaaible (in the caae of constant paraaetera) to derive 
a univorsal cooling law for thla reglae. Thla law has its siaplest fora when t 
it expressed ae a function of T , and an explicit expression for t(T ) is given 
fo ' the caae of arbitrary (but constant) theraal paraaaters (cf. eqs. 26 and 
27). It!) the fora of the predictive problea, the radiance E(T^) may be 
regarded as a forcing function that is producing a "forced decay" in T^(t). 

The "fvee decay modes" arising from an arbitrarily specified interior tempera- 
ture field have died out, and the expression for the interior teaperature 
T(r,t) can be cha .cterized by a fora-lnvarlant functional dependence on the 

forcing fun' iwn E(T ) or, equivalently, on T . That is, 

8 8 

T(r,t) - T[r, T^(t)]. (3) 

S ... an ixpression f it T(r,t) is given, which has the fora of the sum of 
polynomials in r multiplied by powers of T^(t) in which all coefficients are 
tiise-in'.. pendent numbers (cf. eqs. 14 and ?3). This solution corresponds to 
the ''particular solution" in the theory of linear differential equations, but 
the present problea is nonlinear because the forcing function E(T^) is a 
function of the unknown surface temperature*. 



2. Outline of Solution Method ftploy ed 


A four-step procedure is »^>loyed to solve the weak fom of the predictive 

problea: Q) First the sphere is imagined to be iMersed in a controlled 

temperature bath with a given tiae-dependent teaperature T^(t). This is the 

only boundary condition Imposed (except for the requlreaent that the interior 

teaperature be everywhere finite). This is a linear problem whose solution can 

be expressed as a convolution integral. This integral is then transformed so 

as to cast the solution into a slsiple and novel fom In which the temperature 

is expressed as an infinite (but subsequently truncated) sum of characteristic 

polynomials In r of degree (2n+l), each of which is multiplied by (“D^d'^T /dt*'. 

s 

Thus the tine derivatives of the surface teaperature can be regarded as 

"driving" the solution for the interior temperature. (2) the next step consists 

of imposing the boundary condition that relates (3T/or) to E(T ) . When the 

s .s 

solution to the linear problem Is substituted Into this boundary condition, the 
boundary condition assumes the form of an ordinary differential equation In 
T^(t) that Is linear except for the fact that the forcing function E(T^) is a 
nonlinear function of the dependent variable The differential equation Is 

of infinite order, but Is truncated to fourth order. (3) An Iterative pro- 
cedure Is used to find a self-consistent solution to this nonlinear differential 
equation In which dT^/dt and the higher derivatives are expressed as sums of 
powers of This Is the "particular solution" that characterizes the post- 

transient regime. (4) The equation dT /dt - F(T ) is Integrated to give an 

s s 

explicit expression for t(T ), which is the universal cooling curve for a hot 

s 

sphere in the post-transient regime. 

The saae ordinary differential equation that results from the imposition 

of the botmdary condition involving the radiance E also serves as the basis of 

the solution of the calorimetric problea, except that E(t) and T (t) 
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(and hanca all of ita darlvatlvaa In tha ordinary dlffarantlal aquation) auat 
now bo ragardad as known functions, so that tha dlffarantlal aquation bacoaas 
an algabralc aquaclon In which tha unknowns ara a and c. If o and c ara conatants, 
than It la only nacassary to know tha valuas of E and at two aoaanta during the 
cooling in order to construct and solve a syatsa of two algabralc aquations for 
tha two unknowns a and c. Howavar, an altamatlva procedure for solving for 
a and c la employed that leads to axprasslons for than that land thansalvas to 
an easily automated Iterative solution. 

Tha problem of finding the temperature dependence of variable parameters 
o(T) and c(T) Is solved by introducing new dependent and Independent variables 
that reduce the nonlinear partial differential equation to tha same form as 
that for constant parameters except for the presence of a very small effective 
heat source density whose effect can be neglected. The sasM solution procedure 
is followed as In the case of constant parameters, and the end result Is 
expressions for o(T) and c(T) In terms of E(t) and T^(t) (and their derivatives 
that can be solved by a simple Iterative procedure. 

There appears to be no discussion of the weak form of the predictive 
problem In the existing literature, emept for the limiting case In which tha 
Interior temperature may be considered to be constant throughout tha sample. The 
strong form of the predictive problem, l.e. tha Initial valua heat problem, 
has of course a superabundant literature. ^ 

The calorimetric problem has been Investigated by a small group of 
mathematicians^' , but the Investigations have been lladted either to two 
constant parameters, or else only one variable parameter, and the resulting 
expressions are complicated and difficult to solve nuswrlcally. One Important 
cause for this complication Is the fact that these authors Insisted on 







maintaining tha raquiramant of uniform initial interior temperature, and the 
concept of the poat-tranaient regime waa not introduced. Without limiting 
the analyaia to thia regime, and aimultaneoualy rigidly maintaining specified 
initial conditions, the problem becomes overdetermined, and in general ill-posed. 

3. Predictive Problem for Conetant Thermal Parameter s 

Although the primary purpose of this report is to present the solution of 
the calorimetric problem, the necess-'-ry inalysls must be confronted already in 
the simplest predictive problem (l.e. the one for constant parameters). The 
solution of this problem also provides estimates of the numerical values of 
various quantities that become observational data In the calorimetric problem. 

The partial differential equation governing the flow of heat In a body 
In the absence of any Interior heat sources is 

V • (k 7 T) - c 3T/3t (4) 

where k Is the thermal conductivity and c is the heat capacity per unit volume. 

(Strictly speaking, since c Is the product of the mass density and the constant- 

pressure specific heat, and because the mass density Is a function of temperature, 

c Is the heat capacity of the quantity of material that occupies unit volui^* 

at some given reference temperature , and is not strictly the heat capacity per 

unit volume at an arbitrary temperature, although this simplified nomenclature 

will be used In what follows. If p(t) Is the mass density, then the constant- 

pressure specific heat C ■ c/p. C and p enter Into the formalism only as 

P P 

the product c ■ P^p-) 

The boundary condition that must be applied to eq. 4 Is obtained by noting 
that the Integral of the right side of eq. 4 over the volume occupied by the 
hot sphere Is equal to the rate of decrease of heat content of the sphere, and 
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•o Bust bs squsL to 4irR^E whsrs ths rsdianes E Is ths rsdlstsd host flux psr 
unit arss. Applying Gauss' Thaoran to tba voIubm Integral of tha laft slda of 
aq. 4, ths following boundary condition rasults: 

-k(3T/9r)^ - E(t) (calorimetric problaa) (5a) 

■ eoT^ (predictive problem) (5b) 

where the subscript s designates quantities evaluated at tha surface, o Is the 
Stefan-BoltSBann constant, and e Is the emlttance. This boxmdary condition 
asstanes that the sphere is radiating Into a "cold vacuum", l.e. a vacuum at 
absolute sero of teisperature. This assumption Is convenient, but In no way 
essential for what follows. 

With proper scaling, all quantities can be made dimensionless, and all 
parameters (Including dimensionless ones) can be eliminated from the equations, 
reducing them to a universal form that applies to spheres of all sites with all 
possible (constant) thermal parameters. In order to do this. It Is necessary 
to introduce a characteristic length r^, time Inter /al ti, and temperature T^. 

The first two are self '‘■evident : 

ri = R ; tj = R^/a. (6a, b) 

The definition of the characteristic temperature T^ la not so obvious. Various 
definitions can be found In the literature, such as the highest or lowest 
temperatures encountered In the problem, or an equipment-related tesq>erature 
such as that of a furnace to which the sample Is exposed. If, however, one 
wishes to reoiove even dimensionless constants from the equations, only the 
following definition will suffice: 

Ti = (k/Rea)l/3. (6c) 

This characteristic temperature Is determined by the slse and thermal properties 
of the sample, and Is Independent of the sti^rtlng conditions of the experlswnt 
or of the properties of the environment to which the ssmple Is exposed. For a 
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tungsten sphere for \^ich R * 1 cm, ■ l'’,300 K which Is nearly five times 
the melting temperature (3653 K) , whereas for a sphere of uranium dioxide of 
the same size, Tx * 3650 K, which is only 500 K higher than the melting tempera- 
ture (3151 K). 

The dimensionless varla* les (designated by a circumflex) corresponding to 
the characteristic magnitudes deilned above are 

r = r/R ; t = at/R^ ; f = T/Ti ; C/a.b.c) 

and the corresponding equation and boundary condition can be written as follows: 


3^(rT)/3f2 _ 3(^T)/3t - 0 


(8a) 

- (3T/3r)^ - E(t) 

(calorimetric problem) 

(8b) 

- T 

8 

(predictive problem) 

(8c) 

where 



E(t) = E(t) (k(Ti/R)]"l - RE/kTi. 


(9) 


The first expression for E above shows that It may be regarded as a dimensionless 
heat flux for which the scaling flux is [k(Ti/R)], l.e. the flux driven by the 
temperature gradient ti/R. 

As noted In Sectl:.n 2» the first step In the solution of eq, 8 above 
consists of replacing the boundary condition (8c) by the requirement that at 
the surface of the sphere f be equal to the (as yet unspecified) temperature 
T^(t) which Is regarded as the forcing function that drives the Interior 
temperature T(r,t). It Is more convenient to deal with the tentperature 
difference (T - and to this end the dependent variable is defined as 

0(r,t) = ?lf(?,t) - tg(t)]. (10) 
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Thus the equations to be solved are 


(lla) 

(llb) 


320/3^2 - 30/3^ • i dT^/dt 

0(0, t) ■ 0(1, t) ■ 0 . 

This is a linear inhonogeneous equation with Dirlchlet-type boundary conditions 
that is driven by the forcing function r dT^/dt. It is easily confirsned that 
the solution can be written in the following form: 

CO 

0(r,t) X sin(mvf) A (t) exp(-(mit)2t] (12a) 

n VI Ifl 

n»*l 


where 


A (t; » 
n 


t 



expt(mit)2fl [dT^(t')/d€’] dt'. 


(12b) 


An alternative form results from splitting the integral of (12b) Into two parts 
as follows: 

A„(t) - + A (t) (13a) 

ID ID ID 


where Is a ccnstant given by 


o 

r 


*»■ j 


exp[(arir)^t' ] [df^(t')/dt'i dt* 


(13b) 


and 


'■'“■I 


e*p[(miT2) i'] [df (t')/dt'l dt* 


(13c) 


If eq. 13a is substituted into eq. 12a, we have the more familiar form of the 
solution in which the terms involving B describe the rapidly decaying normal 
modes, and the summation involving A^(£) describes the inhomogeneous part of 
the solution that is driven by the forcing function r df^/dt on the right side 
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of eq. 11a. Because A(0y ■■ 0, It is evident that onJy the part of the solution 
(12a) that Involves the B Is non-zero at t ■ 0, and that I(-l)™ B /m] are 

ID ID 

the Fourier coetilclents for a sinusoidal expansion of B(r,0). The longest-lived 
normal mode (m»l) decays like exp (-ir^t) ■ exp l-(7r^a/R^)t] , so when R / x, 
this longest-lived mode has a magnitude of order exp(-TT^) » 5 x 10*^, so it is 
safe to say that by this time the post-transient regime has become established. 

When the solution (12a) is written using (12b) rather than (1?), the 
entire solution is regarded as “driven**, with thi effect"? of the driving force 
reaching back to t ■ - «>. That is, one effectively takes t = - » to mark the 
beginning of the experiment, and takes the position that for all finite £ the 
decaying homogeneous solution has long since died away, leaving only the 
inhomogeneous driven solution. For a solid sphere wnos'^ interior temperatures 
can be influenced only by the environment to which its surface is exposed, this 
is a physically more natural point of view than that posed by the initial-value 
problem, which requires that the initial interior temperature field be 
specified, rather than the past history cf the surface temperature - 

If the derivative in the integrand of (12b) is expanded in powers of 
(t* - t) , the integration can be carried out explicitly and the solution 
assumes the form 

00 

0(f,t) - (-1)" T (t) P„{^) (14a) 

s n 

n»l 


wht ;e 


and 


(t) E d” T /dt" 

8 

rrt 

Pn(r) 5 2 sinCimT?), 


(14b) 


(14c) 
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The are polynomials of degree (2nM) that have the following properties: 


P (0) • 0 for n > 0; 

(15) 

n — 


P„U) - 1. 

(16a) 

P (1) » 0 for n > 0; 
n 

(16b) 

d2p /dr2 - 0, 
o 

for 0< r < 1 

d^P /dr2 ■ - P , for n > 0 . 
n n-*i 

\ U7b) 


The first few of these polynomials can be found in handbooks. (Sec. 
example, p. 87 of reference 11.) Property (17b) can then be used to generate 
successively higher-order polynomials by double Integration using properties 
(15) and (16b) to determine the two constants of Integration. The four 
lowest-crder polynomials are listed below: 


P^,(^) - i (18a) 
Pl(^) - ^(l-r2)/6 (18L) 
Pa(r) - J (7 - 10 r2 + 3 r**)/360 (18c) 
P 3 (r) - r(31 - 49r2 + 21 r** - 3 r®)/15,’20 (18d) 
P 4 (r) - f(381-620 f2+ 29A r** - 60 r® + 5 r®) /1, 814,400. (18e) 


Using the property (17b), it is a very easy malcer to show by direct 
substitution that the differential aquation (Da) is satisfied. The properties 
(13) and (16) guarantee that the boundary conditions (lib) are satisfied* 

The form of the solfition given in (14) represents a different way of 
stating thac the interior temperature field is determined by the history of 
of the surface temperature. Whereas the form of the solution given in (12) 
represents this history by an integral reaching into the distant past (with 
a kernel that provides “memory loss”) , the form given in (14) represents this 
history by meons of all of the time derivatives of the surface temperature 
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calculated at the present Instant of time. Instead of the ever-diminishing 

influence of the surface temperature as one reaches deeper into the 

past, it is now the ever-diminishing signific nee of the hlgher-crder 

time derivatives of surface temperature as the order of the derivative 

increases. This decreasing influence of as n increasns results 

from the tact that P decreases as n increases. For this reason, a 

n 

truncation of the infinite series in (14a) is feasible in the case of a 

nume ’.cal calculation. In the calculation described in this report, 

only the terms of (14a) up to and including n * 4 were retained* 

The infinite set of polynomials does rot span the entire class 

of continuous functions on tho interval £ 1* The set is, however, 

sufficiently large to include all permissible solutions to (11). To see 

this, it should be noted that since (6/f) = f - f , the set (P /r) 

s n 

ought to be capable of describing the most general form cf f - 

inside the sphere, and on the surface should vanish. This latter requirement 

is satisfied in view of (16b) and the fact that n = 0 is excluded from 

the summation in (14a). As to the completeness of the et it 

should be noted thci: (as (18) indicates) it is isomorphic to the set of 

all even powers of r, which corresponds to the requirement that V t = 

0 at the center of the sphere, which in turn follows from the fact that, 
in the absence of a point source or sink at the origin, the heat flux 
must vanish there. 

The second step of the solution method that was outlined in Section 
2 above consists of substitutirg (14) into the boundary condition (5b) 
in order to obtain an ordinary differential equation that determines the 
time dependence of T^(t). The dimensionless fonr of the boundary condition 
(5) is 
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-O0/3r)^j^ “ E(t) (caloriu9!trlc problem) (19a) 

“ T\ (predictive problea) (19b) 

Differentiation with respect to r of the solution given in (lA) , and 
subsequent substitution into (19) yields the follow ing result : 

oo 

^ ^ (-l)*^ I2^°B^/(2n)I] d° T /dt° ■ E(t) (calorimetric problem) (20a) 
n«l 

> f (predictive problem) (20b) 

s 

where are the Bernoulli numbers, vh:.ch are tabulated in sundry handbooks^^. 
The first 7 of these are 


Bi « 1/6, B2 - 1/30, P.3 - 1/42, B4 » 1/30, B5 - 5/66, Be - 691/2730, 
By - 7/6. 

If the series in (20) is truncated at n * 4, the follo%rt.ng ordinary 
differential equation results: 


, d**! 

1 8 

4725 dF*“ 


_X_ ■*'*8 
945 dP 


, d2f , df 

, 1 s 1 s 

45 dF~ “ 3 df 


E(r) (calorimetric problem) 
f ** . (predictive problem) 


( 21 ) 


(22a) 


(22b) 


In the calorimetric problem, T^**^(t) and E(t) are given functions of 

tine, so when (22a) is written in dimensional form it becomes an algebraic 

equation in the two unknowns a and c. This case Is discussed in Section 

4 below. In the predictive problem however, one must solve the ordinary 

differential equation (22b), which would be linear except for the fact 

that the forcing function is a nonlinear function of the dependent 

variable f . Except for this nonlinearity, one could assert that the 
s 

solution of (22b) consists of the sum of a homogeneous part that requires 
a knowledge of four iLiitlal conditions, and a particular solution that 
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is entirely deterained . by the forcing function and is independent of the 

initial conditions. The fora of (14a) shows that Che initial conditions 

of (22b) determine the initial interior teaperature field, and vice 

versa. This hoaogeneous solution has completely died away by the time 

the systea has entered the post-transient regime that is completely 

determined by Che forcing function f**. The "particular solution" chat 

s 

describes this regioe can be found analytically by an iterative process 

which leads to a set of self-consistent expressions for the derivatives 

as fxmctions of f alone, i.e., there are no arbitrary degrees of 
s s 

freedos that must be specified by initial conditions or any other condition. 
These self-consistent expressions are the following: 


t( 1). -[3 f** - ^ T^ + ^ 
s I ® 5 s 25 


ilo.m 

s 875 


f(2)= 36 ^ TlO + ^ t13 , 

s s 5 s 25 s ’ 


t(3)^ -(756 T^O- T^^ , 

S L ® ^ ®J 



(23a) 

(23b) 

(23c) 


f (*"5-22680 T^3 , (23d) 

s s 


If these expressions for substituted into (14a), 0(r,t) 

becomes the sxm of products of powers of r and in which the coefficients 

are mmhers , rather than functions of £• That is, 9 becomes a form- 

invariant function of r and Tg (or, equivalently, of the forcing function 

f , and its total time dependence is contained in the time dependent of 
s 

T (t). This is one way of characterizing the post-transient regime, 
s 

The expressions given in (23) are valid only if the truncation of 
the infinite series in (20b) is valid. The range of validity of this 
truncation is indicated by Table I below which uses the eqs. 23 to 
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calculate the magnitudes of the various terms of eq. 22b for a range of 
values of t\ (Powers of 10 are Indicated In parentheses following the 
numbers, e.g. 3.9(-2) - 3.9*10“^.) 


Table I 


T 

8 

T** 

8 

(_tO)/3) 

(4T<2)/45) 

8 

(-2tO)/945) 

(4T<‘*>/4725) 

1 

1.000 

0.525 

0.384 

-4.709 

4.600 

3/4 

0.316 

0.223 

3.96(-2) 

-6.00(-2) 

0.114 

1/2 

6.25(-2) 

5. 64 (-2) 

4.70(-3) 

+7.92(-4) 

5.86(-4) 

1/4 

3. 91 (-3) 

3.86(-3) 

4.72(-5) 

+1.43(-6) 

7. 15 (-8) 

1/8 

2. 44 (-4) 

2.44(-4) 

3.80(-7) 

+1.48(-9) 

8.73(-12) 


A left-to-rlght sweep of the ^Irst row of Table I shows that the magnitudes of 
the higher-order terms In (22b) are deflnlte'*y not small compared with the 
magnitude of the forcing function f**. However, this Is the case for T ^ 1/2, 
so It seems safe to conclude that In this rSnge the truncation of the Infinite 
series Is a valid procedure. 

It is hoped that It will be possible tj extend the free-cooling calorimetric 
method up to temperatures that exceed the melting temperatures of refractory 
materials by several hundred degrees Kelvin. If the methods end results of 
this report (and Its Intended extension to Include melting) are to be used. 

It Is necessary that the size of the sample sphere be chosen so that the 
condition 

f = T /Ti - T (k/Rej)"^^^ < 1/2 (24) 

8 8 8 “ 

Is satisfied. Referring to the remarks made In Section 1 above. It is evident 
that this condition Is easily satisfied for a tungsten sphere of radius 1 cm. 

In fact, at the melting point of tungsten (3653 K) , ^ 0.211 for a sphere 
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of this size. However, for s sphere of uranlua dioxide of the sst.e size 
at its aeltiag point (3151 K), * 0.86. However, for a sphere of 

radius 1 on, at the aelting point of uraniua dioxide > 0.40, so 
condition (24) is satisfied. Thus, in applying the ftee-coollng method, 
one of the first design considerations is to choose a sample sphere size 
that is small enough to satisfy condition (24). In general, for metals 
eny radius up to 1 cm would be safe, but for refractory non-metals, 1 mm 
should be regarded as the upper limit. 

The eq. 23a can be written in the form 


dt - - dT /F(T ) 
8 8 


(25) 


where F^r^) is the expression contained in the brackets on the right 
side of (23a). This equation can be integrated to give t(f^). When 
this is carried out to an accuracy consistent with that of eq. 23a. the 
following equation results: 


i i (i;’ -1) - «n f. 


(26a) 


where 

R 1/3 

t = ot/R2 and ^ T^. (26b, c) 

The arbitrary integration constant In (26a) has been chosen so that t « 

0 for * 1. It should be noted, however, that (26a) Is reliable only 

for £ 1/2. The first term on the right side of (26a) Is the well- 

known expression that characterizes the cooling of a very good conductor. 

The other terms allow for the presence of a temperature gradient within 

the cooling material. The last two terms are not important. If they 

are omitted, the error that is introduced at f • 1/2 is less than 1%, 

8 

and considerably smaller than that for smaller values of . Eq. 26 is 
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• un'lV3rsal cooling curve titet it valid for apharicel aaoplaa of all 
sires tor aati rials having constant thenal paranaters. It is useful 
for astiaating the cooling tiae At ■ t^ - t^ that is required for a 
sample to cool down froa an initial teaperatura to a final taaq>aratur«3 
T^. Dropping the last two tiras in (26a), this becoacf 



For example, using parameters appropriate to a tungsten sphere of radius 
1 cm, one con calculate the time required for the sphere to cool froa 
3000 K to 1500 K. In this case (27) yields At ■ (501.2 0.62) seconds. 

This shows that the second term, that allows for an interior temperature 
gradient Is negligible, although it would have a relatively greater 
significance in the case of a short time Interval near the beginning of 
the cooling process. This cooling time of 502 seconds is to be compared 
with the time R^/a *3.3 seconds that is required for the post-transient 
regime to become established. The same calculation for a small sphere of 
uranlxtm dioxide having a radius of 1 nmt yields At ■ (28.5 + 0.37) ■ 28.9 
seconds for the sphere to cool from 3000 K to 1500 K. This Interval is 
to be compared with R^/a ■■2.0 seconds, which is the approximate time 
required for Che post-transient regime to become established. 

4. Calorimetric Problem for Constant Thermal Parameters 

In Che calorimetric problem, E(t) and T(t) arc given functions, and 
a and r. are unknown constants. A knowledge of the fact that the cooling 
obeys the t '* radiation law is necessary only for a calculation of the 
eiLictance c, which is given by 

- E(t)/o lT(t)]‘» , ■ (28) 
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where a Is the Ste£sn>9oltXMUu coastant. This relstioc Is actusUy the 
deflolclon of e» sod it Is valid even If e Is not constant. Fron a 
knowledge of e(t) and T(c), one could find the teaperature* dependent 
ealttance e(T). 

The basic equation is now (22a) rather than (22b). When the deflnltione 
contained in (7) are used to convert (22a) into dlaenefonal form, it can 
be written as follom: 



- T (1-F) - 3E/RC 
s 

(29a) 

wh^re 




T^ s d T^/dt 

(29b) 

and 

' 15aT ^ 

(29c) 


This yields the follotd.ng expression for c: 

c - [3E/R (-T^)l (1 - ¥r\ 0:>) 

For a perfect conductor, F • 0 because a ■ «* . For finite a, F -► 0 
as the sphere cools down because the cooling curve T^(t), although still 
descending, tends to flatten out so that Che higher derivatives of T^(c) 

Chat appear in F tend to vanish. This suggests a slnq>le procedure for 
evaluating c and a: First, corresponding readings of E and near the 

cold end of the cooling curve are used to calculate c from (30) with F * 

0. Then this value of c together with corresponding readings of E and 
near the beginning of the post-transient regime are used in (29a) to 
calculate a numerical value for F. Using this numerical value of F in 
(29c) , together with the known numerical values of the derivatives of 

T (t), a cubic equation in (1/a) results which can be solved for a. 

® • 

This value of a could then be used in (30) together with the cold tesq^erature 
readings to calculate an l^iroved value for c, and so on. 
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A slaplc analytical axpraaalon for o can bo derived by neglecting 
and In (29c). To thia approzlawtlen, 


P - 


ISoT 


•• »2 j • 

lt‘" <-V- 


(31) 


and when thia is used in (29a), the following expreaaion for a results: 

3E 


^ . djstk'^ 

15q L dt 


(32) 


RcT 


This expression for a involves c on the right side, but this can be 
eliminated by a substitution. If (29a) is dlffereatiatad, and 
and are Ignored as they were In (31), as well as T^, the following 
approximate expression results: 

(33) 


- - 3 E/Rc. 


With this, (32) becoBies 


15a L dt J ■ L dl J 

It can be shown that the following exrxt expression can be derived 
without making any of the above approxlauitlons: 

k ■ - * <»’v] - - «<“’«] • 


(34) 


(35a) 


where 


f(a,T ) H 
s 


1 

10 


ii i-llM 

a T **‘3 ’ 


8 

d^E 

2 r2 1 

s 

d^E 

149 

R*» 

1 

d‘*l 

b 

IF “ 

45 a E 

IF ■ 

472r 


E 

dt»J 


T = dT /dt. 

8 8 


(35c) 

(35d) 
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Because. this expreseion involves only the ratios T and It 

is Insensitive to the calibration of T and E, which la a significant 
experlaental advantage, lurthemore, because this expression Involves 
only a and not c, ar iterative process can be eaq>loyed to calculate an 
accurate value of a using readings of E(t) and T^(t) and their derivatives 
at a single point on the cooling curve. Then these sane readings and 
the calculated value of a can be used to calculate c from (30) without 
the need for further iteration. In fact, the F that appears In (30) can 
be written in the following form: 


F - 


15a 


dfn (-T ) 
s 

dt 


- f (o 


• v] 


(36) 


where f(a,T^) Is defined In (35b). 

The factor (1-F) that appears In (30) Involves derivatives of T . 

8 

it Is possible to Slake use of (35a) to exchange this factor for one that 
Involves only derivatives of E. The resulting expression for c Is 


where 


c - [3E/R(-Tg)](l+G), 


(37a) 

(37h) 


where g(a,E) Is defined in (35c). It seems likely that, as a matter of 
experimental technique, T^ ana Its derivatives will have higher precision 
than E and its derivatives, so (30) Is likely to give better accuracy 
than (37). 

AHhough these expressions for a and c appear to Involve readings 
taker, at only one point on the cooling curve. In actual fact they Involve 
readings at several neighboring points, because these are necessary to 

calculate the various derivatives of T and E which are not measured 

s 
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directly in the experlaenti However, once theee derlvetivee heve been 
celculeted (end any neceeeery eaoothing of the date carried out) , the 
expreeeions (35) and (30) (or (37)) can be used to march along the 
cooling curve and calculate an Independent pair of valuea for o and c at 
each point of the curve. Thua, if theee thermal parametera are actually 
mild functione of temperature rather than constants, this procedure 
ought to give reasonably accurate expressions for a(T) and c(T). 

There are two sources of error in this procedure, however, which 
will be corrected in the next section. One is that the calculated 
values of a and c should not be associated with the surface temperature 
of the point on the cooling curve at which they were calculated, but 
rather with a temperature resulting from some kind of averaging throughout 
the volume of the sphere. Second, the temperature dependence of the 
thermal parameters makes the basic differential equation nonlinear, and 
cot Just the boundary condition. In spite of these considerations, the 
application of the results of this section will provide good first 
approximations for Insertion into the formulas that are derived in the 
next section. 

5. Calorimetric Problem for Temperature-Dependent Parameters 

The strategy that will be employed in chls section is to invoke 
appropriate redefinitions of both the dependent variables in order to 
reduce the problem to one which, in its dimensionless form, is identical 
to the dimensionless form of the calorimetric problem for constant 
parameters. In this way a formal solution of the equation and boundary 
conditions is achieved. When the transition is made back to the dimensional 
quantities appropriate to the experiment, the effects of the changes in 
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variable* chat vara aada at tha baglimiiig of cha solution aanlfast 
tbaaaalvas, vltb tba result that tanw Involving Che derivatives of the 
paraaatars vltb respect to teaperature appear In the expressions for a 
and c. These expressions are the desired generalisations of those 
derived In the preceding section. This approach autonatlcally refers 
everything to the surface teaperature T^ and to the values of the paraaeters 
at the surface, even though the effects of the higher teaperatures of 
the Interior are lapllctly taken Into account. 

The equation and boundary condition are the saae as given In eqs. 4 
and 5a, except that now k(T) and c(T) are arbitrary and Independent 
functions of teaperature. Instead of working with T is the dependent 
variable, the enthalpy density h(T) now assuaes this role where 


h(T) 



c(T') dT’ + h^. 


(38a) 


and 

hr H h(T^). (38b) 

fs any convenient constant reference temperature. Strictly speaking, 
h(T) Is actually the enthalpy of that quantity of matter which occupies 
unit volume at the reference temperature, but which at a different 
temperature occupies a slightly different volume because of the temperature 
dependence of the mass density. For simpllctly, however, h will be 
called simply the enthalpy density. In terms of h, the partial differential 
equation and the boundary condition become 



V • Ia(T) V h] - 3h/3t 

(39a) 

and 


- Oh/3r)^.j^ - E(t) 

(39b) 

where 


a(T) = k(T)/c(T) 

(40a) 
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and 


a = a (T ). 

a 'a' 


(40b) 


a la a function of t via the tine depcndancc of T (t) , and thla function 
8 8 


o^(t) will be used to define a dlaenalonless time variable t: 


T = R”2 J' a^(t') dt’. 


(41) 


Using (41), the eq. 39a can be put Into the following form: 
V2 h + V • V hj - R~23h/3f 


(42) 


where the second term on the left Is an effective heat source density 
arising froa the effective heat flux -[(a/a^) - 1] This effective 

flux vanishes at the surface of the sphere where a ■ a^. Because of 
this, the volume Integral of the effective heat source le zero. If the 
effect of this effective source density were taken Into account by an 
Iterative perturbation technique. It would be found that h would have to 
be corrected by a small amount 6h such that ^6h would be very nearly 
equal to the effective flux density ■*[(a/a^) - 1]^. Because this vanishes 
at the surface, so would V6h and (36h/3r)^. Because of this, the correction 
5h would produce no change In the boundary condition (39b) which gives 
rl»c to the equation from which the expressions for a and c are derived. 

In other words, the second term on the left side of (42) would produce a 
slight change In the Interior h-fleld, but It would produce no significant 
effect on the equation that Is derived from the boundary condition 
Involving E, and this equation Is all that matters for the solution of 
the calorimetric problem. 

The reference temperature T^, and the values of the sundry quantities 
at this temperature, are used to make the equation and boundary condition 
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dlitens ionless. Dlneaslonless quantities will be designated by a circumflex. 
Thus 

c - c/c^ = c(T)/c(T^); and a = o/ot^ = o(T)/o(T^). (43 a,b) 


The characteristic temperature Is now defined In terms of the values of 


the material constants and at the reference temperature: 
Ti= (k^/Re^o)^^^ 

This Is used to scale both h and T: 


h = h/c 




c(T') dT' + hj./c^Ti 


(44) 


(45a) 


where 


T = T/Ti. 


The dimensionless Independent variables are 

t 


r a r/R and 


T = o R~^ I a (t') dt* . 

' J • 


(4Sb) 


(46a, b) 


In conq>lete analogy with the variable 0 defined in (10), it is convenient 
to introduce the variable n defined by 

n (r,x) = r |ii(r,T) - h^(r)j . 

In terms of these variables, the equation becomes 

S^n/af^ + f [y - d hg/dxj - 3n/3T 


(47) 


(48a) 


where 


--'^1 

art 8 3r J 

3f L“s arj 


(48b) 

(48c) 
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The boundary conditions are 


- On/3i)j,^ - E/a^ (50) 

where 

“ 9 = a(Tg)/a(Tr); ® (51a, b.c) 

Except for the presence of the ^ term In (48a) , this equation and the 
homogeneous boundary conditions (49) are Identical to the system defined by 
( 10 ) and (11), and could be solved In the same way if the iy term were Ignored. 
The Inclxislon of this term presents no fundamental problem since it can be 
handled by an Iterative approach that treats it as a prescribed forcing function 
like (-dh^/dr), except that y depends on r as well as t. The formal solution 
of (48a) Including y(^>'c) (regarded as a specified function rather than as a 
function of n) has been carried out, but will not be glve.a here, because 
numerical analysis (not yet carried out) will almost certainly show that the 
error caused by Ignoring the ry term in the case of any material of practical 
Interest Is not greater than the error caused by truncating the Infinite series 
foL n. For this reason, the ry term will be ignored, and the solution for 

A 

n(r,r) will be taken to have the same form as that given In (14). 

The same argument that converted the boundary condition (19a) Into the 
ordinary dlffsrectlal equation ( 22 a) now leado to the following equation: 

, d**fi , d^h , d^fi , dh 

1 S 2 8 1 8 1 8 /co\ 

4725 dr"* 945 dr^ 45 dt^ 3 dt * ^ ^ 

This la Identical in form to (22a) except that T^, t, and E are replaced by 
h^, t, and E/a^ respectively. To convert to the dimensional form tt is only 
necessary to note that It follows from (45a) and (41) that 
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dh dT /,, 

di " dT \d1 

Successive dlfferei 
taken to note that 

dc 


c(T ) 

B^c 

dT 

8 

^ s 

8 

c Ti 
r ^ 

ot c Ti 

dt 

s r ^ 


in the 

same way. 

with 


(53) 


dT 




dT 

dt 


2 dc dT 

dT dt 
8 8 


(54) 


and similarly for da /dr. 

s 

The same procedures employed in Section 4 that led to expressions 
for a and c in the case of constant parameters aov lead to similar expressions 
that Include the effects of the variability of the parameters. The expression 


for a H a (T ) is 
s s 


15a 


TdT-- 




) - ^ 


d£nE , 

dT- « 


3" 


where 




dJ.no 
8 

dt 

2 _ ^ 
21 o 


8 

d2jn 


dJn c 
s 

dt 

‘ d^T 
1 s 

• dt^ 
T 

s 
c 


10 a 


s T 


, d^T 
1 s 

• dt ^ 


s 




d^JZ.n a 

s 

~d^ 
d£n c 


din a \ 


din a 
£ 

dt 


dt / 


and 

s 


(55a) 


1 / 

^ dt^ 

8 

/dJn 0 

Hit" - ~/V 

“s . 1_ ^ (1 ^ 1 d“E 

dt 35 o E dt? ■ A5 o E dP , 

8 I .8 / \ 

*■ '’“.no filn a \ 

dt ^ \ dt / J 


d^Jn o_ 
dt" 


d Jn E 
- 3 dt 


(55c) 
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where the teres involving eC*) and third and higher derivatives of 


in o and in c as well as cables in the first derivatives of these have 
8 8 

been neglected. 


The two alternative expressions for c = c(T ) are 

3£ (, r 2 *. * ,ir‘ 

'• R(-T ) ( L 


(56) 


where given by (55b), and 


3E 


R(-T^) 


'1 + 


15a 




(57) 


where g' (E,a .) is given by (55c) . 


The point of view to be taken in applying these expressions is that 

E(t) and T (t) are giver, and after suitable smoothing their sundry time 
s 

derivatives are calculated. These functions are then substituted into 

(35) and (30) to obtain the tine-dependent functions a(t) and c(t) which 

are taen used for a and c on the right sides of (55) and (56) to 
s s 

evaluate nore accurate functions a (t) and c (t) which are reinserted 

s s 

into the right sides of (55) and (56) , and so on until the process 
converges. From the final functions cig(t) and Cg(t) , together with 
T (t) , the functions r»(T ) and c(T ) can be found. Because a and c are 
functions of 1 alone, and the functional relations are the same in the 
ia^'eriot of the sphere as on the surface, the subscript s can be dropped 
since the computed functions are the desired functions c:(T) and c(T). 

It is probable that a numerical analysis will show that, for the 
kinds of functional relations a(T) and c(T) that are encountered in the 
case of any real material of practical interest, many of the terms in 
(55b) and (55c) are negligible in the siuise that they are smaller than 
the experimental error of the larger terms. 
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6 . 


Discussion of Fart B 


The working formilas for the application of the free-cooling calorimetric 
method are (55) and (56) (or eq. 57). It is noteworthy that (55) (the 
expressicn for a) is completely insensitive to the units used to measure 
E and since only relative values enter into the formula. Thus, for 
the calculation of a, it is not necessary to calibrate the Instruments 
in absolute terms since arbitrary units can be used. In the case of the 
calculation of c, th^ value of (E/f^) is the only quantity that must be 
calibrated. If tue value of c is known for t) e sample in question at 
any temperature (most likely the coldest temperature in the experimental 
range), then the calibration can be made at that tenq>eratvre. If c is 
completely unknown for the material in question, then an entirely different 
material and different sphere size can be used for calibration purposes 
if c for the calibration material is known at some temperature. 

The derivation of the working formulas made no assumption idiatever 
about t^'"* dependence of E on T^. Rather, it was simply assumed that 
E(t) and T^(t) were independent given functions. Only in the calculation 
of the emittance ^(^g) was any assumption made about the physical nature of 
the cooling. (Cf. eq. 1.) Thus the formulas for a and c apply also to 
a sphere that is cooled by some mechanism other than radiation into a cold 
vacuum so long as the assumption of spherical symmetry remains valid. 

Theti. is no fundamental reason why the methods enq>loyed in this report 
could not be adapted to cylindrical geometry. This geometry admits the 
possibility of controlled electrical heating of refractory materials (an 
established experimental technique). Such heating would be represented in 
the formalism as a prescribed heat source density. The definition of the 
post-transient regime would have to be modified only to the extent of 
including the prescribed heat source as one of the determining factors. 
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along with the cooling law and the thenal paraaetera. Such prescribed heat 
sources have already been Included In the analysis for spherical geometry. 
Although their Inclusion conplicates the fomallsa. It Introduces no 
fundaaental difficulties. 

Following a numerical test and confirmation of the relations derived 
In this report, the proposed next stage of development of the analysis Is 
the Inclusion of the possibility of a solid spherical shell surrounding a 
molten core that shrinks In size as the sphere cools. The object oz this 
analysis would be to determine the heat of fusion, l.e. the Jump In enthalpy 
density across the solidification front. For this purpose h Is a much more 
natural variable than T, so the formalism of Section S should provide a natural 


starting point. 
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Task 5 


Experimental and Theoretical Studies in Wetting and Multilayer Adsorption. 

M. R. Holdover and J. W. Schmidt 
Thermophysics Division 
Center for Chemical Physics 
J. W. Cahn 

Center for Materials Science 
Summary 

Equipment has been assembled and/or fabricated to measure the 

ellipsometric parameters of the liquid-vapor interface of binary mixtures 

near the predicted transition from high adsorption to low adsorption. 

Preliminary measurements on the mixture C-,F,. - C-H, , have not yet discovered 

7 14 7 14 

the predicted transition. A variety of interfering phenomena have been 
tracked down. Further measurements are in progress. 
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Introduction 


tfe have assembled a seMi-autonatic ellipsoneter together with 
saaple cells and a thermostat capable of precisely monitoring the ellip- 
sometric parameters of the liquid*vapor interface above binary liquid 
mixtures. Ve have made a computer model of the ellipsometer in order to 
optimize its sensitivity to the appearance of the predicted high ad- 
sorption layer. We have tested a method of measuring and recording the 
index of refraction of the binary mixture under study simultaneously 
with the recording of the ellipsometric parameters. The index of re- 
fraction measurement is extremely sensitive to the appearance of macro- 
scopic amounts of a second liquid phase; thus it will serve to locate 
the adsorption transition with respect to the miscibility gap in the 
specific liquid sample under study. We expect the strategy of simul- 
taneously measuring the index of refraction and the ellipsometric para- 
meters to obviate the necessity of carrying out detailed exacting 
chemical analyses of each sample of each mixture to be studied. 

1 2 

In preliminary measurements we have searched for the predicted * 

high adsorption-low adsorption transition in the liquid-vapor interface 

between a C_F, . - C-H. . mixture and air saturated with this mixture *s 
7 14 7 14 

vapor. We have found no evidence of such a transition in a sample which 
undergoes phase separation near 23^0. (The consolute temperature in 
this system is near 46^C.) This sample is approximately 7% fluorocarbon 
on a mole fraction basis. In these preliminary measurements, the ellip- 
sometric parameters showed a smooth temperature dependence between 23^ C 
and 60^ C which would be equivalent to the adsorption of a 70A thick 
layer of the fluorocarbon at the interface between the vapor and this 
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hydrocarbon-rich liquid. This toaporaCure dependance is probably 
apurious. At the aoaant our highest priority is to detemine its 
origin. These preliainary aeasureaents could have detected an abrupt 
high adsorption-low adsorption transition equivalent to a 20A change in 
the thickness of a fluorocarbon-rich layer at the liquid-vapor inter- 
face. We expect laproved sensitivity in future measurements and to scan 
the comrosltion range in order to locate and trace out the expected 
interface transition curve on a phase diagram. 

The index of refraction measurement revealed a secular drift in the 
temperature at which bulk phase separation occurred (roughly 0.5°C/mouth). 
This was traced to a leak in the quartz sample cell which we have now 
repaired. 

The measuresient system we are developing will be automated. This 
is required because of the long equilibration times encountered. 

Experimental Procedures 

The experimental procedures we have carried out will be described 
under the following subheadings: (1) Optical Table, (2) Thermostat, (3) 

Equilibration (4) Sample cell, (5) Model for ellipsometer, (6) Optical 
conponents for ellipsometer, (7) Optical components for xndex of re- 
fraction measurement, and (8) Sample preparation. When it is appropriate, 
we will indicate the rationale for the procedure followed as well as the 
improvements we anticipate in the near fucure. Figure 1 shows the 
experimental arrangement. 

(1) Optical Table: The elllpsometric measurement is sensitive to 

the angle of the incident light. As the liquid-vapor surface sloshes 
about (from ambient vibrations) the angle of incidence varies. Further- 
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■or*, th* path of tha raflactad baaa through various optical conponants 
is also a function of tha instantanaous position of cha surface. These 
two effects ■odulata tha intensity of the raflactad bean. Thus they 
sppaar as noise in tha allipsoaatry. A coaBercially nanufactured air 
suspension table was purchased and Installed to use as a low vibration 
support for the entire expcrisient. Even at the conparltively low 
sloshing frequencies the table significantly reduced the liquid-vapor 
surface vibrations. The tsble is helu very nearly level by a systen of 
valves which adait coapressed air to the legs. When functioning pro- 
perly the vsriations in level are + 1 ailliradian. This seeas to be 
satisfactory at the present precision of our measurements. We have 
found that occasionally one of the valves may stick. Then, much larger 
excursions in the position of the table level occur. These excursions 
are correlated with room temperature variations (which alter the air 
pressure in the table suspension) and are corrected periodically as the 
valve finally opens. This larger amplituc'e tilting produced an un- 
acceptably large ellipsometric signal whose origin was tricky to locate. 
Ultimately we nay have to regulate the table level with a much smaller 
dead band than the coupressed air servo system tolerates. This is 
possible by moving two small weights on the table. 

(2) Thermostat: We are now using a two stage thermostat. The 

temperature of the outer stage is controlled by circulating water from 
an external circulater /thermostat. We have used this circulator between 
IS C and 60°C. The temperature of a particular point on the outer shell 
is now monitored with a thermistor. Temperature fluctuations of this 
point are on the order of + 0.05®C. 
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The therul relaxation tlae of the aample la several hours (see 
below). We change the teaperature setpoint of the outer stage by switch- 
ing reoljtors In parallel with the teaperature sensitive element of the 
regulator. The switching is done at preselected times. This hardware 
based prograaaing scheme will eventually be replaced with a software 
based one incorporating feedback. 

The inner shell of the thermostat 5s regulated by a home-made 

electronic servo system using a thermistor sensor. The set point 

-4 

temperature exhibits noise of less than 10 C; however, there are sub- 
stantial temperature gradients across the inner stage (roughly 30 mK) • 

In the near future two improvements will be made to reduce this gradient 
by nearly two orders of magnitude. First, the point heaters on the 
inner shell will be replaced by a distributed heater placed symmetrically 
around the shell. Second, the temperature of the outer shell's set 
point will be servoed in such a way as to maintain the outer shell's 
temperature much closer to that of the inner shell. 

(3) Equilibration: Equilibration of concentration gradients in 

binary mixtures via diffusion is a notoriously slow process. Typical 

-5 2 

diffusion constants in the liquids of Interest are 10 cm /s. 

Typical dimensions of our samples are 3 cm; thus characteristic times 
should be several hours when one liquid phase is present and much longer 
when two liquid phases are present. A magnetic stirring bar has been 
enclosed in a glass envelope and placed within the sample cel . This 
bar is turned periodically (e.g. 2 minutes every 30 minutes) by a second 
magnet built into the thermostat which in turn is attached to a motor 
outside the thermostat. We are able to directly observe equilibration 
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of the ellipeoBetrlc paraaetere of the llquid*vapor surface a«.ter a 
temperature change. When one liquid phase is present, equilibration 
takes about one hour. When two phases are present, equilibration can be 
monitored using the index of refraction measurement because che Index of 
only the upper phase Is measured. We have observed that equilibration 
took between two and four hours when the sample was cooled into the two 
phase region even though the liquid-vapor interface acts as a nucleation 
site for the lower fluorocarbon-rich liquid phase and the sample Is 
stirred vlgourously. 

(4) Sample cells: The sample cells are essentially hollow prisms. 

They are fabricated from fused quartz windows and quartz tubing held 
together with a high melting temperature solder glass. The cells are 
filled through a quartz tube wh ch Is sealed with a teflon stopper and 
0-rlng. (It Is possible to seal the cells by fusing the tube shut; 
however, at this point we have not done so). The prism angle Is 60^ and 
Is oriented with its apex upward. The angle of Incidence for the 
elllpsometer Is near 60". This minimizes the Influence of the residual 
birefringence of the cell’s windows on the elllpsometer performance. 

This angle of Incidence is also close to Brewster’s angle for the 
hydrocarbon-rich mixture; thus. It Is near optimum for the elllpso- 
roeter’s sensitivity. 

In equilibrium, the liquids under study are In contact with quartz 
and vapor (the stopper Is in the vapor phase). The hydrocarbon-rich 
phase of mixtures totally wets the quartz. Thus the sample 

cells' walls and corners as well as the glass encapsulated stirring bar 
do not act as nucleation sites for the appearance of the fluorocarbon 
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rich phAse. The fluorocarbon-rlch phase can only appear at the liquid- 
vapor Interface, the Intel face we are studying elllpsonetrlcally. The 
absense of "traps*' for the fluorocarbon-rlch phase Insures that the 
Index of lefractlon measureoient will clearly Indicate the initial 
formation of macroscopic amounts of this phase. 

(5) Model for the Elllpsometer: The ellipsometer arrangement Is a 

conventional one: source-polar Izer-compensator-sample-analyzer-detector. 

In order to choose reasonable declgn compromises, a computer program was 
written to calculate the intensity at the detector as a function of the 
angle of Incidence, the orientations of the polarizer, compensator, and 
analyzer, the index of refraction of the substrate fluid, and the 

thickness and index of refraction of the multilayer adsorbed film we are 

searching for. This program Is an essential tool for the design of the 

instrument and the interpretation of the measurements* With it we can 

explore the sensitivity of the measurements to extraneous factors such 

as the tilt of the table and component imperfection. 

(6) Optical components for the elllpsometer: A ImW He-Ne laser is 

used as a source. The polarizer and analyzer crystals have an extinction 
ratio near 10^. A Faraday modulator is Inserted between the polarizer 
and analyzer. In effect, it rotates the plane of polarization through a 
few tenths of a degree. The modulator is excited with an AC signal so 
that the photomult iplyer output can be monitored with a lock-in amplifier. 
A DC signal is also fed back to the modulator to insure that the Intensity 
of the detector is near null. The ellipsometer is set up for extinction 
when the sample is at a high temperature. Then we expect essentially 
zero preferential adsorption at the liquid-vapor surface. As the sample 
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is cooled the DC conpcnent of the current fed back to the modulator Is 
continuously monitored. If a thin dielectric film, l.e. an adsorbed 
layer, appears on the liquid surface the extinction will be maintained 
by the feedback signal. This method does fail when the layer becomes on 
the order of 1/8 of a wavelength of light thick. Then two optical 
components must be adjusted to maintain extinction. 

An integrator with a 2 minute period is used in the feedback loop 
to average over the interface vlbraL .ons. 

The most serious problem we have encountered with this arrangement 
Is an ambient temperature dependent residual elllptXclty Introduced by 
the modulator. We have started monitoring this by periodically inserting 
a test sample between the modulator and a second detector. This sample 
simply a quartz flat whose optical properties are conr'-ant. In the 
future, we will attempt to use thermostated modulator of our own design. 

(7) Optical components for index of refraction measurement: The 

index of refraction measurement is based on measuring the deviation that 
a collimated beam of light experiences upon passing through the prism- 
shaped sample cell. A second He-Ne laser Is u ;ed as a light source. 

Its beam Is spatially filtered and expanded to fully Illuminate the 
cell. The beam passes through the cell and Is reflected back upon 
Itself. (This doubles the sensitivity of the Index of refraction measure- 
ment and reduces the number of windows needed for the thermostat). A 
beam splitter Is used to separate the returning beam from the Incident 
one. A telescope focuses the returning beam upon a position sensitive 
detector. The voltage from the detector provides a continuous monitor 
of changes in the Index of refraction of the sample. 
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(8) Sample preparation: The mixtures ar*? prepared from 

the pure componerts received* A rough measurement indicates the 
consolute temperature of our samples ars within PC of that reported in 
the literature for these components. Because ve intend to establish the 
location of the adsorption transition in relation to the phase separation 
temperature it is not necessary to prepare mixtures of a precisely kr wn 
composition. Instead we prepare mixtures such that two phases are 
present at the temperature at which phase separation Is desired. Then 
the hydrocarbon-rich upper phase is drawn off and used as tne sample. 
Experimental Results 

As we have indicated, the experimental techniqujs are still being 
actively developed. The tests of the apparatus that have been carried 
out to date have not yielded new results of physical significant. We 
have not yet observed the predicted phase transition, but we have 
scanned only one composition. 

Con c lusions and Discussion 

The measurement technique we have developed can be improved sig- 
nificantly. Substantial progress is expected in the near future. An 
automatic aata acquisition system has been ordered. As of this writing, 
the primary factor limiting the resolution of our present equipment is a 
large temperature dependent shift in the apparent etllipsometric para- 
meters of our fluid samples. 
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